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Abstract. The Frequency Domain Reflectometry (FDR) is studied as a powerful 
tool to detect hard or soft faults in star-shaped networks of nonuniform lossless trans- 
mission lines. Processing the FDR measurements leads to solve an inverse scattering 
problem for a Schrodinger operator on a star-shaped graph. Throughout this paper, 
we restrict ourselves to the case of minimal experimental setup corresponding to only 
one diagnostic port plug. First, by studying the asymptotic behavior of the reflection 
coefficient in the high-frequency limit, we prove the identifiability of the geometry of 
this star-shaped graph: the number of edges and their lengths. The proof being rather 
constructive, it provides a method to detect the hard faults in the network. Next, we 
study the potential identification problem by inverse scattering, noting that the poten- 
tials represent the inhomogeneities due to the soft faults in the network wirings. Here, 
the main result states that the measurement of two reflection coefficients, associated 
to two different sets of boundary conditions at the extremities of the tree, determines 
uniquely the potentials; it is a generalization of the theorem of the two boundary spec- 
tra on an interval [3] . 
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1. Introduction 

The rather extensive literature concerning the "inverse scattering problem" and the 
"inverse Sturm-Liouville problem" on graphs have mostly followed separate pathways 
except for a very few results [3, [6] . In the following paragraphs, we briefly recall the 
previous results on these subjects and at the end we will situate the result of this paper 
with respect to the others. Indeed, as it will be seen later, the inverse Sturm-Liouville 
problem considered in this paper raises from the necessity of finding a minimal setup 
for solving the inverse scattering problem. 

The paper [8] considers a star-shape graph consisting of N infinite branches and 
solves the inverse scattering problem assuming the measurement of N — 1 reflection 
coefficients. Next, in the paper [9], Harmer provides an extension of the previous 
result with general self-adjoint boundary conditions at the central node. This however 
necessitates the knowledge of N reflection coefficients. 

The paper [13] studies the relation between the scattering data and the topology of 
the graph. They show that the knowledge of the scattering matrix is not enough to 
determine uniquely the topological structure of a generic graph. In pQ, Avdonin and 
Kurasov consider again a star-shape graph with N finite branches. They, show that the 
knowledge of one diagonal element of the response operator allows one to reconstruct 
the potential on the edge corresponding to this element. 

As mentioned above, in parallel to the research on inverse scattering problems, 
another class of papers consider the inverse problem for Sturm-Liouville operators on 
compact graphs. These results can be seen as extensions of the classical result provided 
by Borg [3], on the recovering of the Sturm-Liouville operator from two spectra on a finite 
interval. The main progress in this field has been made by Yurko [IE] , [Hj and [20] . The 
paper [18] deals with the inverse spectral problem on a tree. The idea is to generalize the 
Borg's result in the following sense: for a tree with n boundary vertices, it is sufficient 
to know n spectra, corresponding to n different settings for boundary conditions at the 
extremities, to retrieve the potentials on the tree. In the recent paper [20], the same 
kind of result is proposed for a star-shape graph including a loop joined to the central 
node. Finally, [19J provides a generalization of [18] to higher order differential operators 
on a star-shape graph. 

In [15], the author proves that under some restrictive assumptions on the spectrum of 
a Sturm-Liouville operator on a star-shape graph with some fixed boundary conditions, 
the knowledge of this spectra can determine uniquely the Sturm-Liouille operator. 

In this paper, we consider the inverse scattering problem motivated by the 
application in fault-detection/diagnostic of star-shape LC transmission networks. We 
are interested in minimal experimental setting providing enough information on the 
network and the potentials on its branches. The graph consists of N finite branches 
joined at a central node and we add a infinite branch to this central node for the 
experimentation. We will see that the knowledge of the reflection coefficient is equivalent 
to the knowledge of the spectra for the Sturm-Liouville operator defined on the compact 
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part of the graph for various boundary conditions at the central node. By considering the 
associated inverse Sturm-Liouville problem, we will show that under some assumptions 
on the geometry of the metric graph, the knowledge of only two reflection coefficients, 
corresponding to two settings for the boundary conditions at the terminal nodes, is 
enough to determine the potentials (at least locally). 

In the next section, we explain the application under study and the associated 
experimental setup. We will show that the inverse scattering problem for the so- 
called Telegrapher's equation (lossless case) on the network is equivalent to an inverse 
scattering problem for a Schrodinger operator over the metric graph of the network. In 
Section [3l we will consider the direct scattering problem and we will characterize the 
reflection coefficient in terms of the fundamental solutions for Sturm-Liouville operators 
on branches. In Section HJ we will show that the knowledge of only one reflection 
coefficient is enough to identify the lengthes of the branches of the metric graph. This 
result will be useful to locate hard faults (open or short circuits). Finally, in Section [5] 
we consider the main problem of recovering the potentials from the knowledge of one or 
two reflection coefficients. This result will be useful to locate soft faults (local variations 
of the electrical characteristics). We will prove the equivalence of the inverse scattering 
problem with an associated inverse Sturm-Liouville problem on the compact part of 
the graph. This inverse problem is then treated applying the methods extending the 
classical result by Borg [3]. 

2. Frequency domain reflectometry 

The electric signal transmission through a wired network is, generally, modeled with 
the "Telegrapher's equation" and characterized by the parameters L, C, R, G (functions 
of the space position z along the transmission lines) representing the inductance L, 
capacitance C, resistance R and loss conductance G per unit of length. These parameters 
allow a rather complete and understandable description of the transmission lines and 
are sufficient to represent the lines in the frequency range used during reflectometry. In 
the sequel we will suppose that this model can be used for all frequencies. However, 
it appears to be impossible to retrieve all these parameters uniquely through the 
information provided by reflectometry experiments. Everywhere, through this paper, we 
will consider the simpler nonuniform lossless situation (R = G = 0). As we will see later, 
the reflectometry measurement is still not enough to retrieve the both parameters L and 
C but rather an aggregate of these two parameters, the local characteristic impedance 
Z c {z) := y/L(z)/C(z). 

Following [12] and [10], the presentation of the reflectometry experiment of 
this section, is based on a model derived from the "Telegrapher's equation" and 
parameterized by Z c . To cope with the network case, we have translated the Kirchhoff 
rules at the nodes of the network within this new modeling framework. Note that, in 
this paper, we restrict ourselves to the case of a simple star-shape network and therefore 
the only node of the graph where the Kirchhoff rules need to be adapted is the central 
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one. The faults, in which we are interested through this approach, are represented by 
the lengths of the branches (hard faults) and by the heterogeneities of Z c along the 
branches (soft faults). The considered reflectometry experiment model is based on a 
far-field method consisting in adding a uniform infinite wire joined to the network at its 
central node. In practice, connecting a matched charge to the extremity of a finite line, 
is sufficient to emulate the electrical propagation through an infinite line. 

The linearity of the transmission line model allows to replace any test by an 
equivalent test in harmonic regime. We can therefore start by stating the Telegrapher's 
equation in the harmonic regime, i.e. the tension and the current intensity are 
respectively of the form e~ tuJt V(u, z) and e~ tuJt I(u), z), where u is the time frequency 
and z the position. On each line, we have 

d 

—V(u), z) — iujL(z)I(uj, z) = 0, 

a (1) 

— I(u, z) - iuoC{z)V{uo, z) = 0, 

Through this paper, we assume that 

Al the distributed parameters C(z) and L(z) are twice continuously different iable on 
the transmission lines; 

A2 they are strictly positive, C(z) > 0,L(z) > 0; 



A3 the characteristic impedance Z c (z) := yL{z)/C{z) is continuous at the central 
node of the star-shape network; 

A4 the transmission lines are uniform in a neighborhood of the extremities of the 
branches. 

The Liouville transformation. Note that the reflectometry experiment leads to 
observing the tensions and currents along the time at some position: only the travelling 
times (and amplitudes) of waves are accessible by such experiments. A fault can only 
be localized in terms of the traveling time of the reflected test wave starting from the 
test point. This leads to a particular change of variables, the Liouville transformation, 
allowing to work with the traveling time rather than spatial coordinates. Let us recall 
this transformation: 

x(z) = [ y / L(z)C(z)ds 
Jo 

which corresponds to the wave traveling time from the position to the position z. 
Remark that after this transformation, u is also the wave number on each branch. 

The inverse transformation being well defined, we will write C{x) = C(z(x)), 
L(x) = L(z(x)), V(u,x) = V(u),z(x)) and I(u,x) = I(u,z(x)). 

The Telegrapher's equation (TjQ) becomes 

d 

— — V(u, x) — iuZ c (x)I(u>, x) = 0, 

d , (2) 

— x) — wjZ c (x) 1 V(lu, x) = 0. 

(y JU 
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The wave decomposition and equivalent forms of the Telegrapher's equation. Define 

1 v ( ^_ x dZ c {x) 

dx 

following decomposition of V: 



V±{uj,x) = \{V{w,x) ± Z c (x)I(u,x)) and Q{x) = -Z c {x)- 1 £ ; . We have the 



V(u, x) = V + (uj, x) + VL.(w, x), 

d _ (3) 

— V±(u, x) =F vu)V±(u, x) = ±Q(x)(V + (uj, x) - V-(u, x)) 

In particular, in an interval where a branch is uniform, Q(x) =0, and the solution 
is the sum of waves of opposite directions. For any x and x in this interval: 

V(u,x) = V + (to,x)e tuj{x - s) + V„((j,x)e- tuj{x - s) . 

Define now y(u>,x) = ^Y c (x)V(cu, x) and q(x) = y/Zjx)-^ y/YJjc) , with Y c {x) = 
Z c (x)~ l . The Telegrapher's equation (CQ) becomes a Schrodinger equation: 

d 2 y 

- -j^(uj,x) +q(x)y(uj,x) = u 2 y(u,x). (4) 

It can be seen from (jHJ) and (jlj), that, the knowledge of the potential Q(x) or of q(x) 
and of the boundary conditions on V and I is sufficient to compute the solution of (PQ) on 
the network. In our lossless situation, we have chosen q as the parameter to be identified 
through the reflectometry experiment. A variant of ([3]), the Zakharov-Shabat equations, 

would be the good choice in the more general lossy case. Remark that — Q 2 + q = 0. 

dx 

The reflection coefficient. Taking e tu ( x ~ 1i \ oriented toward the increasing x, as 
a reference forward wave, the reflection coefficient is the following ratio of backward 
over forward wave amplitudes: R(u,x) = v~(ux)e-*» x = e 2 ^ x ' ^~|^ . In particular, 
R(u,x) = R(u,x): defined in this way, the reflection coefficient is constant in intervals 
where Z c is constant. For an arbitrary x, denoting by Z(u,x) = V(u,x)/I(u,x) the 
(possibly infinite) apparent impedance at x, we still define the reflection coefficient as 
being 

RM=e^ T X) - Z 7 fy (5) 
Z(uj, x) + Z c (x) 

With this definition, in general, V-(u,x) = e~ 2tuJX R(u,x)V + (uj,x), and, using Q, it is 
easy to check that R is solution of the following Riccati equation: 

9B &> X ) - e - 2wx Q(x)R(LU, x) 2 + Q{x)e 2wx = 0. (6) 
ox 

Finally, note that, we will only consider "positive real" terminal impedances Z(u,t), in 
the sense that Z(— u, r) = Z(u,r) and %tZ(u), r) > 0. This together with the fact that 
Z c [x) G M. is positive implies that \R(uj,x)\ < 1. Furthermore, experiments with uj > 
are sufficient. 

The network under test. Throughout this paper, F represent the compact star- 
shape network consisting of the branches {ej)^ =l joining at the central node and r + is 
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the extended graph where the test branch e is also added to the graph. We have N + 1 
equations of the form 

d 2 y 

- + Qj(x)yj = u?yj x G (0, Tj), (7) 

where Tj is the wave traveling time associated to the branch number j (to = oo 
as the added branch eo is assumed to be an infinite line). In particular note that, 
as the infinite branch eo is assumed to be a uniform transmission line, we have 
qo(x) =0, x E (0, oo). 

The boundary condition for the reflectometer. Consider now that a generator with 
a matched internal impedance Z c and an electromotive force 2V g e tuJX is connected in x of 
an interval X where Z c is constant. We have V(lu,x) + Z c (x)I(lu,x) = 2V g e tulx which can 
be the boundary condition if the branch terminates at x. We have V + (u,x) = V g e lu)X 
and V(u,x) = V g {e tuJX + R{u,x)e~ tuJX ) for all x e 1 . As, R(u,x) = R(u,x), the 
reflection coefficient can then be determined from the measurement of the tension 
V(uj,x) anywhere in X. In the sequel, we will use a test branch eo with a constant 
Z c connecting a matched generator to the central node . The measured reflection 
coefficient on this branch, will be simply written R(u). As we will work with y, we 
choose V g = a/ Z c (x), so that y(uj,x) = e luJX + R(uo)e~ tulx . Finally, it will be convenient 
to take the same positive orientation on all the branches, from the central node at x = 
toward the increasing x. Our reference forward wave on e is then in the direction of 
the decreasing x, so that, changing x into —x, and supposing e of infinite length, the 
boundary condition for the reflectometer is: 

y(u, x) ~ e~ luJX + R(u)e wx as x -> oo in e . (8) 

The two sets of boundary conditions at the network extremities. In order to recover 
the potential of the star-shape network, we will need to consider two experimental 
settings, with open circuit or short circuit at the extremities of the branches. This will 
lead to a problem similar to solving an inverse spectral problem for the Sturm-Liouville 
operator when two spectra are known. 

The first setting corresponds to open circuit configuration at the extremities of the 
finite branches ((e-,)^). This, together with the Assumption A4 on the local uniformity 
of the lines around r/s, leads to boundary conditions of the form Ij(cu,Tj) = 0, or 
equivalently, we obtain the setting called, the Neumann configuration: 

2/;.(r,) = j = l,...,N. (9) 

The second setting corresponds to the short circuit configuration at the extremities 
of the finite branches ((e.,)^). This leads to boundary conditions of the form 
Vj(u,Tj) = 0, or equivalently, we obtain the setting called, the Dirichlet configuration: 

y j (T j ) = J = 1,---,AT (10) 

Remark 1. In some of the applications that we have in mind, the reflectometry 
experiment has to take place without perturbing significantly the normal utilization 
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of the transmission network, so that using open or short circuits conditions would be 
impossible. There is a way to circumvent this problem by computing the results of the 
open or short circuit experiments from results of two less invasive experiments. The 
idea is to use nonlinear superposition properties of solutions of Riccati equations as in 
[17] . in order to get a closed-form representation of the reflection coefficient, solution of 
(jHJ), as a function of a general load impedance (value of Z at the extremity of a branch) 
and of two particular solutions corresponding to two load impedances more compatible 
with the network utilization. This will be presented in a forthcoming paper. 

The boundary condition at the central node. It writes 

yi(u,0) = Uj(u,0) =: y(u) i, j = 0, • • • , N, 

2^2/i(w,0)-2/ (w,0) = -- z0 

j=l c 

where y'j (u,0) and (Z^')'(O) denote the spatial derivatives at the point x = and Z 3 C is 
the characteristic impedance of the branch number j. Note, in particular, that we have 
applied the continuity of Z^s at the central node (Assumption A3): Z 3 c {jS) = Z®, Wj. 

Formulation of the model. In conclusion, in order to study the LC-transmission 
line equations on the graph T + , we can study the Schrodinger operators 

d 2 

D(C^ V ) = closure of C% v in H 2 (T + ), (12) 

where C^(T + ) (resp. C^(T + )) denotes the space of infinitely different iable functions 
/ = ®^=ofj defined on Y + satisfying the boundary conditions 

./•«') ./••(()) j,/ = 0,---,A, 



E ft® - /o(0) = Hf (0), H = -- V zQ L (13) 

j=l c 

/j(ri) = (/ J (r J )=0for^(r + )), j = !,•••, A. 



3. Direct scattering problem 

The operators (C^ v , D(Cjj- v )) are essentially self-adjoint. To prove this fact we observe 
first that these operators are a compact perturbation of the operators ®" =0 ( — j with 
the same boundary conditions. Now, we apply a general result by Carlson [4J on the self- 
adjointness of differential operators on graphs. Indeed, following the Theorem 3.4 of pE], 
we only need to show that at a node connecting m edges, we have m linearly independent 
linear boundary conditions. At the terminal nodes of {ej}^ =l this is trivially the 
there is one branch and one boundary condition (Dirichlet or Neumann). At the central 
node it is not hard to verify that (llip define iV + 1 linearly independent boundary 
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conditions as well. This implies that the operators (C^- v , D(Cj^ v )) are essentially self- 
adjoint and therefore that they admit a unique self-adjoint extension on L 2 (T + ). 
We are interested in the scattering solution where a signal of frequency uj is applied at 
the infinite extremity of the infinite branch. In such a case, we will be seeking a solution 
satisfying the asymptotic behavior 

y (x, uj) ~ e~ luJX + R(uj)e l " x , for x -> oo. 

Proposition 1. Under the assumptions Al through A4, there exists a unique solution 

continuous with respect to uj, of the scattering problem, satisfying 

• ~~£zyj( x i u ) + 1j( x )yj( x ^ w ) = uPVjfa for j = 0, ■ ■ • , N; 

• (yj(x,uj))j =0 satisfy the boundary conditions (Tl\) and (G|) or ( T70j) ; 

• For each uj G R, there exist R(uj) such that 

y {x, uj) ~ e~ luJX + R{uj)e iujx , x -> oo. (14) 

VKe wz// denote the reflection coefficient R(uj) defined by (fl^[ ) m i/ie Neumann (resp. 
Dirichlet) case by Rj\f(uj) (resp. Rd{uj)). This coefficient appears to be unique. 

Proof. This proof gives us a concrete method for obtaining scattering solutions. Indeed, 
we will propose a solution and we will show that it is the unique one. 
In this aim, we need to use Dirichlet /Neumann fundamental solutions of a Sturm- 
Liouville boundary problem. 

Definition 1. Consider the potentials q 3 as before and extend them by on (— oo,0) 
so that they are defined on the intervals (—00,7^]. The Dirichlet (resp. Neumann) 
fundamental solution (p J v (x,uj) (resp. <p J j\f(x, uj)), is a solution of the equation, 

d 2 



rVW x ' <*>) + q^x^Jx, uj) = uj f^Jx, uj), x G 



"OO, T; 



^ v {Tj,Uj) = 0, (^)' v (Tj,Uj) = 1, 

vm^j,^) = i, <y )V( r i> u ) = °- 

Consider, now, the function 



where 



$! Vj jv(x,oj) = ® j=o y J vM (x,u), 
^° vM (x, uj) = e~ tuJX + R vM (uj)e wx , x G [0, 00), 
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Here the coefficients Rv,jv and o? v ^ are given by the boundary conditions (1111) at 
the central node: 

l + R Vf M-(u) = c4) >Af {uj)(fi)^(0,Lo), j = l,---,N, (15) 

N 

Y J ^ vM {u){i ( P v ^)'{Q,oj) + - ife^M) = H(1 + Rdst(u)). (16) 

3=1 

One, easily, sees that this ^v,Af is in D(£j^ v ), the domain of the operator, and satisfies 

the conditions of the proposition as soon as the coefficients Rv,M and {pSpjJ) are 
continuous with respect to u>. This, trivially, provides the existence of a scattering 
solution. Here, we show that ^t>,jv is actually the unique one. 

Assume that there exists another Y-djj = ®jLo^r>,Ar( x i u ) solution of the scattering 
problem. By the compact injection of H 2 in C 1 , we now that Yp^x, u) and ^/ VJ ^(x, u) 
are C 1 functions of x over [0, 7}]. Therefore, the Wronskian 

is well-defined. Moreover, as the potentials qj(x) are continuous functions over [0, Tj] and 
as Ypfj-(.,uj) and ^ 3 t> k[{-, are solutions of the associated Sturm-Liouville equation, 
they are in fact C 2 over [0, Tj]. Thus, the derivative of the Wronskian is also well defined 
over [0, Tj]. Through a simple computation and by noting that Yp u>) and ^ 3 T) u-{.,uj) 
are solutions of the same Sturm-Liouville equations, one has 

-^W(Yi M (x, u), %^{x, u)) = 0, x G [0, Tj ], 

and so the Wronskian remains constant over the interval [0, Tj]. 

For the finite branches {ej)^ =1 , applying the (Dirichlet or Neumann) boundary 
conditions at the terminal nodes, we easily have 

W{Yl^{ T j,uj),¥ ViM {Tj,u)) = Q, 

and therefore the Wronskian is identically over the whole branch. This implies that, 
Y*(.,uj) and % ^(., 10) are co-linear: 

Y vM (x, u) = (3 J v ^{uj)(f vM (x, u), x e [0, Tj], j = 1, • • • , N. 

Over the branch eo, as Ypj^(.,uj) satisfies a homogenous Sturm-Liouville equation 
{Qo = 0); it necessarily admits the following form 

What remains to be shown is that one necessarily has Rv,Af{u) = Rv,Af{^) and similarly 

Indeed, the equations (fT5j) and ffl6l) provide N + 1 linear relations for the N + 1 
unknown coefficients Rv,m and Trivially, as soon as, the coefficients 
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{ ( Pv at(^> u; ))j-i are non " zero ! these linear relations are independent and there exists 
a unique solution for the unknowns R-pjV and («^, A ^)^_ 1 . However, the zeros of each 
one of the coefficients {tfP v ^(0, correspond to isolated values of uj (square-root of 

the eigenvalues of the operator — +qj(x) with Dirichlet boundary condition at x = 
and Dirichlet or Neumann boundary condition at x = tj). Therefore, the coefficients 
Rv,m and (o4?,Af)^-i are well-defined except for a set of isolated values of uj. However, 
as these coefficients need to be continuous with respect to u), they will, also, be defined 
uniquely over these singular points. 

Furthermore, dividing (|T5|) by {1 + R v ^{uj)) and inserting (1T5|) . we find the explicit 
formula 

1 - RvMu) _H 1 A (^)'(O, 
Finally, inserting the value of Rd^{uS) found in (|T7|) into (fT5|) . we find 



□ 



4. Inverse scattering for geometry identification 

As a first inverse problem, we consider the inversion of the geometry of the network. In 
fact, we will prove the well-posedness of the inverse problem of finding the number of 
branches N and the lengthes (Tj)^ =1 of a star-shape graph through only one reflection 
coefficient Rj^(uj) (the case of Dirichlet reflection coefficient can be treated similarly). 
Furthermore, as we will see through the proof of the Theorem [TJ the method is 
rather constructive and one can think of an algorithm to identify the lengthes, at 
least approximately. The proof is based on an asymptotic analysis in high-frequency 
regime of the reflection coefficient and some classical results from the theory of almost 
periodic functions (in Bohr sense). Before, announcing the main Theorem, we need 
a few lemmas. A first lemma precises the high frequency behavior of the Neumann 
fundamental solutions (yr^-)^: 

Lemma 1. Consider a potential q in C°((— oo, r]) and take the Neumann fundamental 
solution, (pj\f(x,u), defined as in Definition^ We have 

<pj\f(0,uj) = cos(cjt) + £>(-), as uj — > oo, 

w 

(<pj\f)'(0, uj) = u sin(a;r) + O (1), as uj — > oo, 
where (<f_\f)'(0, uj) denotes the spatial derivative 4-cpj\f(x,uj) at x = 0. 
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Proof. We start by writing <fjy(x, u) in its integral representation. Indeed, the 
fundamental solution (pj^(x, uj) is given by [H]: 

(prt(x, uj) = cos(lu(t — x)) + / Gj^{t — x, r — s) cos(u;(t — s))ds, 

J x 

where G_\f(x,y) is a real function with the same regularity as jj q(s)ds: here, it is C 1 
with respect to both coordinates. We note that, as Gj\f is in C 1 , by integrating by parts, 
one has: 

r 1 s=t 

/ Gj^(t, t — s) cos(u;(r — s))ds = GV(t, t — s) sin(u;(r — s)) 

Jo u 3=0 

+- f T T-( G ^(^ T - s )) sin (^ - s))ds = o(-\ 
uj J ds \ujJ 

Therefore ipj\f(Q,uj) = cos(ujt) + O (1/co) as u> — > oo and we have the first relation. 
For the spatial derivative (<pf\f) r (x,uj) at the point x = 0, we have: 

((pjf)'(0, uj) = uj sin(c<j(r))— G(r, r) cos(u;r)+ / — [Gj^ij— x, r— s)) cos(u;(r— s))ds. 



dx x=o 



The kernel y) being C 1 , we have 



r d 

— G(t, t) cos{ujt) + / — (Gjs[{t — x, t — s)) cos(u(t — s))ds — O (1) as UJ —>■ oo. 
Jo x=o 

Thus, (</Pj^)'(0, a;) = u;sin(a;(r — x)) + O (1) as uj — > oo and the second relation 
follows. □ 

As we see by Lemma UJ in the high-frequency regime, the fundamental solutions 
become asymptotically independent of the potential q. The next lemma, provides an 
explicit method to identify the number N and the lengthes {tj)^ =1 of the branches for 
the homogenous case, where we know that the potentials {qj)f = i are all zero. 

Lemma 2. Consider a star-shape network T composed of rij branches of length Tj 
(j = 1, • ■ ■ ,m) all joining at a central node so that the whole number of branches N 
is given by Y^j=i n r Assume the potential q on the network to be (q = 0). Then 
the knowledge of the Neumann reflection coefficient Rj^{uj) determines uniquely the 
parameters {nj)^L 1 and (t^)^. 

Proof. We need to apply the explicit computation of the reflection coefficient provided 
by (fTTl) . The fundamental solutions are given, simply, by ipj^(x,uj) = cos(u;(r 7 - — x)). 
Therefore: 

1 — R\f{uj) 1 rT 1 v-^ ujsvairjuj) 

-H — y nj — 



1 + Rj\[(uj) iuj iu cos(rju;) 
The knowledge of Rm{uj) determines uniquely the signal: 



3=1 
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Assuming, without loss of generality, that the lengthes Tj are ordered increasingly 
7"i < • • • < r m , the first pole of the function f(u>) coincides with 7r/2r m and therefore 
determines r m . Furthermore, 

n m = lim cos(W m )/(u;), 

and therefore one can also determine n m . Now, considering the new signal g{oj) = 
f(u>) — n m tan(co>r m ), one removes the branches of length r m and exactly in the same 
manner, one can determine r m -\ and n m -\. The proof of the lemma follows then by a 
simple induction. □ 

We are now ready to announce the main theorem of this section: 

Theorem 1. Consider a star-shape network T composed of nj branches of length Tj 
(j = l,---,m) all joining at a central node so that the whole number of branches 
N is given by ^2J l = i n j- Assume the potential q on the network to be, simply, C°. 
Then the knowledge of the Neumann reflection coefficient Rj\f(uj) determines uniquely 
the parameters (nj)™^ and (Tj) 1 JL 1 . 

Proof. Assume that, there exists two graph settings (r J -,gj)^ 1 and {fj,qj)^ =l (the 
lengthes Tj are not necessarily different) giving rise to the same Neumann reflection 
coefficients: Rj^(cu) = Rj^{oj). By the explicit formula f|T71) . we have 

V (^)'(0, W ) _ if(^)'(o )W ) 



(,i ^—^ ,J Yn ,.,\ 1,1 j—~< 



This is equivalent to: 

JV / N 

n^(o,w) ^(^y(o^)n^(o^) 

3=1 \k=l tyk 

N I N \ 

n ^(o, U ) ^(^)'(o, u) n ^(o, u) = o. (is) 

3=1 \k=l tyk J 

Defining the function: 

N / N \ N I N \ 

F(oj) := Y\ cos(wf fc ) ^sin(cjr fe ) JJcos(urrz) - JJcos(an-j-) ^ sin(wf fc ) cos(a;f z ) , 

3=1 \k=l l^k ) j=l \k=l lj=k ) 

the asymptotic formulas of Lemma [1] and the ffl8l) imply 

F(uj) = 0{1/uj) ascj^oo. 

However, the function F(u) is a trigonometric polynomial and almost periodic in the 
Bohr's sense [2]. The function F 2 (u) is, also, almost periodic and furthermore, we have 

M(F 2 ) := lim ^ / F 2 {u)du = lim ^ ( I F 2 {u)du + [ F 2 {u)du 

Sl^ooiZ Jq C-^oo S I \Jq Ji 

1 „ ] 



< lim - d + C 2 / — du; ) = 0. 
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This, trivially, implies that F = (one only needs to apply the Parseval's Theorem 
to the generalized fourier series of the function F). However, the relation F(u) = is 
equivalent to 



N N 



i=i i=i 

and therefore, by Lemma [2], the two settings are equivalent and the theorem follows. □ 



5. Inverse Scattering for potential identification 

A second inverse problem, related to the detection of soft faults in the network, can be 
formulated as the identification of the potentials on the branches. Here, we consider 
the case of homogenous perfect wires. We will show that the measurement of the two 
reflection coefficients Rv{oj) and Rjyiu), corresponding, respectively, to a short circuit 
and an open circuit experiment, is enough to identify uniquely the small changes in 
the potential remaining in a certain regularity class. Indeed, we will prove that the 
inverse problem of finding the H 1 ^) potentials in an e L°°(r)-neighborhood of the zero 
potential (homogenous case), is well-posed for e small enough. In this aim, we will need 
an additive assumption on the electrical lengths (rj)j =1 , in order to remove symmetries 
leading to degeneracy problems. However, as it will be discussed later, it seems that 
this assumption can be relaxed. 

In a first result (Theorem [2]), we prove that under some natural assumption on 
the electrical lengths, the knowledge of only one reflection coefficient (here Rj^(u)) is 
sufficient to identify uniquely the values: 

Qj = / Qj{s)ds Vj = 1,---,N. 
Jo 

If anyone of these quantities appear to be different from zero, we know that there has 
been a change of parameters in the corresponding branch. By performing classical 
inverse scattering techniques over this branch [T2l [TO] we can identify its soft faults. 
The question is therefore to identify the soft faults in the network which does not 
change the quantities (qj)jL v This is treated in the two Theorems H] and [3j For both 
these theorems, we will need some more restrictive assumptions on the electrical lengths 
{ T j)jLi- I n Theorem El we will show that the knowledge of only one reflection coefficient 
(here Rjy(u)) is enough to identify the potential when we know that it has not changed 
on the first half of the branches. The Theorem H] provides a well-posedness result for 
the inversion of the potential over the whole graph but necessitates the knowledge of 
both reflection coefficients Rx>{oj) and R^{uj). 

In the sequel, we note that the potential over the infinite branch eo is always given 
and is 0. Indeed, this homogenous line is added only for the reflectometry experiment. 
The following theorem provides a global inversion result concerning the quantities q~j. 

Theorem 2. Consider a star-shaped graph Y and assume that 
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Bl The electrical lengths {rj}^ =1 are not entire multiples of each other, 

j- ■ 

— £ N for i^ j. 

n 

If there exist two potentials q = ®^ =1 qj and q = ®^ =1 qj in H 1 ^) giving rise to the same 
reflection coefficient, Rj^ioj) = Rj\t{uj), one necessarily has: 

qj(s)ds= qj(s)ds j = l,---,N. 
Jo 

Remark 2. Note that the result of the Theorem is also valid if we have Rt>{oj) = i?x>(a>). 

This theorem allows us to identify the situations where the soft fault in the network 
causes a change of the quantities q~j. As explained above, the inverse problem can 
then be considered on separate branches and solved through classical inverse scattering 
techniques. In the two following theorems, we assume that the soft faults in the network 
leave the quantities q~j unchanged and as for the perfect situation we are dealing with 
homogenous lines we will assume that 

B2 qj = qj (s)ds = for j = 1, • • • , N. 

Theorem 3. Consider a star-shaped graph F and assume that 

Bl' For any j,j' G {1, • • • , iV} such that j ^ j' , Tjjry is an algebraic irrational number. 

Then, there exists e > small enough such that, if the potentials q and q belonging to 
and satisfying B2, ||<2i|z,°°(r), ||<?IU°°(r) < e and qj(s) = qj(s) for s e [0, Tj/2], 
j = 1, • • • , N, give rise to the same Neumann reflection coefficients, Rj^(cu) = Rj^{uj), 
then q = q. 

Remark 3. Once again, the result remains valid if we replace the Neumann reflection 
coefficient by the Dirichlet one. 

Remark 4. The assumption Bl' seems very restrictive. However, as it will be seen 
through the proof, the only thing we need is that for any of fractions tj/tj/, there exists 
at most a finite number of co-prime factors (p, q) E N x N, such that the Diophantine 
approximation 

IL _ P 
Tj> q 

holds true. However, this is a classical result of the Borel-Cantelli Lemma that for almost 
all (with respect to Lebesgue measure) positive real a's this Diophantine approximation 
has finite number of solutions. 

Theorem 4. Consider a star-shaped graph F satisfying Bl'. There exists e > 
small enough such that, if the potentials q and q belonging to H 1 (F) and satisfying 
B2 and ||qi|L°°(r), ||?IU°°(r) < e, give rise to the same Neumann and Dirichlet reflection 
coefficients, 

Rj\f(uj) = RjKf(ou) and Rv{w) = Rt>(u), 

then q = q. 



< 
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Remark 5. The Theorem H] is a natural generalization to the case of a graph of the 
Theorem of two boundary spectra on an interval [3] (see, for instance, Theorem 1.4.4, 
Page 24, [5]). In [3j, Borg proved that the knowledge of two spectral data corresponding 
to two boundary conditions, determine uniquely the potential on an interval. Here 
the reflection coefficients Rj\f(u) and Rv(uj) play the role of this spectral data (see the 
Subsection 15. ip . 



5.1. Inverse Sturm- Liouville problem 

Throughout this subsection, we will consider a general star-shaped graph V (of N finite 
branches) and a potential q = ®f =1 qj belonging to H l (T). Furthermore we assume for 
the potential q that, the norm ||g||L°°(r) is sufficiently small. 

The main objective of this subsection is to show that the knowledge of the reflection 
coefficient R/^(uj) for C^- (resp. Rv(u) for £j) is equivalent to the knowledge of 
different positive spectra of Sturm-Liouville operators defined on T with Neumann (resp. 
Dirichlet) boundary conditions at terminal nodes and for various boundary conditions 
at the central node. In fact, defining the function 



[1 + Rjw(w)) ' 
where H is given by ( 1131) . we have the following result. 

Proposition 2. Fix w G I and define the Schrddinger operators Ctfp(uj) on the 
compact graph V as follows: 

d 2 

D{C M>v {u)) = closure of C^ V (T) m H 2 (T), 

where C^?^(r) (resp. C^ V (T)) denotes the space of infinitely differentiate functions 
f — ®f=ifj defined on T satisfying the boundary conditions 

/;(0) //(0) ■■/ j,j' = l,---,N, 

N 

j'=i 

/ifo) =0 (/ J -(r i ) = 0/ O rC^,(T)), J = l,---,N. 

Then we are able to characterize the positive spectrum of Cn^u) as a level set of the 
function h^^{oj): 

Proof. We prove the proposition for the case of Neumann boundary conditions. The 
Dirichlet case can be treated exactly in the same manner. We start by proving the 
inclusion 

o- + (C x (u)) C {e | f G R, MO = M«>)} • 
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Let £ 2 G a + (C^(u)), then there exists \I/ eigenfunction of the operator Cj^{ui) associated 
to £ 2 . In particular, it satisfies 

N 

2*;.(o) = m<*>)*, 

where \I/ is the common value of $ at the central node. 

Now we extend ^ to the extended graph r + , such that is a scattering solution 
for Cjj (see the Proposition [1]). In particular, the function must satisfy, at the 
central node, 

*t(0) = *+(0), j = l,---,N, 

N 

^(^)'(0)-(^ + )'(0)=^+(0). 

i=i 

Noting that \l/ is an eigenfunction of (C_\f(uj), D(£n{uj)), we have 

TV 

M"W(0) - W)'(O) = - W) = (0). (19) 

Now, noting that over the infinite branch admits the following form 

= R M {i)e lix + e~* x x G [0, oo), 

the relation ( fl9l) yields to 

/wHOMO + 1) - <(MO - 1) = #(M0 + 1), 

or equivalently 

This proves the first inclusion. Now, we prove that 

<t + (£a/-M) 3{e 2 |eeR, MO = M")} • 

Let ( G 1 be such that TiaKO = We consider a scattering solution of 

the extended operator (defined by ([ED ) associated to the frequency £ 2 . We, then, 
prove that the restriction of to the compact graph V is an eigenfunction of Cj^{uS) 
associated to the eigenvalue £ 2 . This trivially implies that £ 2 G a + (Cj K f(u)). 

In this aim, we only need to show that this restriction of \& + to T is in the domain 
D{C^{uj)). Indeed, this is equivalent to proving that the boundary condition: 

JV 

5>+)'(0) = /^)*o + (0), (20) 
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5>+)'(0) = H*+(0) + (*+)'(0) = (if + (0). 



Furthermore, 



tt+(0) = Rjsit) + 1 and (*+)'(0) = t^AO ~ 1) 



and so 



JV 



*+(o) = MO«tf(o) = /iH^(o). 



□ 



Rm(0 + 1 

This proves (1201) and finishes the proof of the proposition. 

We have shown that, the knowledge of the reflection coefficient Rtfp(uj) for T + is 
equivalent to the knowledge of the positive part of the spectra for £/v",d f° r ah values h 
at the central node. 

The following proposition provides the characteristic equation permitting to identify 
the eigenvalues of the operator Cj\f jT> (uj): 

Proposition 3. The real A 2 > is an eigenvalue of the operator Ctf t v(w) if and only if 
^at,»(A) = fy/vr,x>M$A/VD(A), 

where 



i=i 



i=0 



(21) 



(p^f£>(x, A) feeing t/ie fundamental solutions on different branches. 



Proof. We give the proof for the Neumann boundary conditions, noting that the 
Dirichlet case can be treated, exactly, in the same manner. Assume A 2 to be a positive 
eigenvalue of C^{u). The associated eigenfunction, y\{x) = has necessarily 

the following form: 

where a/s are real constants and the vector (ax, ■ ■ ■ , ajv) is different from zero. The 
function y x , being in the domain D(Cj\f(u)), it should satisfy the associated boundary 
condition at the central node. This implies that the vector • • • , a^) is in the kernel 
of the matrix: 



M 



Mr(0,A) 





-<^(0,A) 
<^(0, A) -<^(0, A) 
<^(0, A) 



\-Mw)y>V(0, A) + ^(0, A) ^2,(0, A) ViV(0, A) 



\ 



W0, A) / 



Short title 18 

where ^-(0, A) denotes -^<pA x ' A) \ x =o- This means that the determinant det(M) is 
necessarily 0. Developing this determinant we find: 

^(A) = /wM«MA). 

□ 

Corollary 1. Consider two potentials q = <%>f=iqj and q = (8>jLi9j and denote by Cj^- 
and the associated Neumann Schrodinger operators defined on the extended graph 
T + . Assuming that the of the reflection coefficients Rj^{uS) and Rm{oj) are equivalent 
Rj\f(u>) = Rtf(u), we have 

$/tf{w)Vtf{u)) = $jv(w)*jv(w), Wu e R, (22) 

where \l/_/v-, ^JV an d are defined through (j2*T|) for the potentials q and q. 

Proof. By Proposition [21 uj 2 is an eigenvalue of the operator Cm{uj) and £jv(a>). 
Applying the Proposition [31 this means that 

^f(u) = hx(ui)<&tf{u) and ^!{uj) = hj^(u)^j^{Lo). 

As Rjj(uj) = Rj\f(uj), we have hj^(ou) = hj^iu) and thus the above equation yields 
to (El. □ 



The above corollary is also valid when we replace the Neumann by Dirichlet 
boundary conditions. Finally, this corollary yields to the following proposition on the 
difference between the two potentials q and q. 

Proposition 4. Consider two potentials q = ®jLi?j and q = (8>jLi?j and denote by Cjj 
and Ctf, the associated Neumann Schrodinger operators defined on the extended graph 
T + . Assuming that the of the reflection coefficients R_\f(u>) and Rj^{uj) are equivalent 
Rtf{yj) = Rj^iyj), we have 

N m 

X)II^ ' W ^W ' W ) / ^A x WA x ^WA x ^) dx ^^ Vwgk, (23) 

3=1 k+j J ° 

where qj = q~j — qj . 

Proof. For j — 1, • • • , N, we have: 

q(x)<f^ r (x,u)(p^ r (x,u)dx- / qj(x)(ffa(x,w)(fo(x,w)dx = 



JO 

= vA x ' u )^A x > u )\ x x=o - ^A( x ^)A( x ^)\ x x=o 

= ^(0, utfAo, u) ~ ^(0, w)^.(0, w). 
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Here the second line has been obtained from the first one, replacing qj(x)(p 3 J y(x,u) by 
-^cfp^x, to) + to 2 Lpj^(x, to) and integrating by parts. Using fl22|) and the above equation, 
we have: 

N p 

3=1 k+j J ° 

N 

II ¥&(°> ^(°> w ) (^(0, w)<^(0, u) - ^(0, ^)^(0, a;)) = 
j=i k^j 

^^(to)®^) - $ M (to)^^(to) = 0. 

□ 

Before finishing this subsection, note that, once more, the above proposition is also 
valid for the case of Dirichlet boundary conditions and Rt>(uj) = Rx>{to) implies: 

N m 

22Yl ( Pv( Q ^)'fU ^) / qj(x)<p } v (x 1 u)0' 13 (x 1 u)dx = O, \/to e R. (24) 
We are now ready to prove the Theorems [2} [3] and HI 
5.2. Proof of Theorem^ 

We prove the Theorem [2] applying the characteristic equation ( l23|) and the integral 
representation of the functions iffo v (x , uj) . Again, for simplicity sakes, we give the 
proof only for the case of Neumann boundary conditions, noting that the Dirichlet case 
can be done in a similar way. 

Similarly to the Section HI the fundamental solutions (p 3 jy(x,uj) are given by [14] : 



<^(x, to) = cos(to(rj — x)) + / G^Tj — x, Tj — s) cos(to(Tj — s))ds, 

J x 

where Gj\f(x,y) are C 1 with respect to the both coordinates. 

The above equation yields to the following representation for the product functions: 

tptf(x, a;)v?j^(x, to) = cos 2 (to(rj — x)) H — / Kfofa — x, Tj — s) cos(2a;(r 7 - — s))ds (25) 

2 J x 

where Ky\f(x,y) is a Volterra kernel, i.e. Kj(x,y) = if y > x and 
Kfox, y) = 2[Gy\ f (x, x-2y) + G^(x, x - 2y)\ + 

/•x px—2y 

+ / Gy\f(x, s)G J j^(x, s — 2y)ds + / Gj\f(x 7 s)G J j^(x, s + 2y)ds. 

Jly—x J —x 

At this point, we note that, as Kjj-(x,y) is a C 1 function, we have 

i J 3 K j H {jj - x, Tj - s) cos(2tu(r i - s))ds = O f^j , x 6 [0, Tj}. (26) 
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Applying the characteristic equation (1251) and developing the products ipjy(x, uS) 
(p 3 j^{x,ui) by (1251) . and finally using (1261) . we have: 

52 ( TTcos 2 (W fe ) ) / ^(s) cos 2 (u;(r i - s))ds = O ( — ), 
o=i Wi J Jo W 



N 



j=i \fc^i 



e (n c ° s2 (^)j \[ (27) 

In the last passage, we applied the fact that J* 3 qj(s) cos2(c<j(Tj — s))ds = O (1/uj), since 
q is in ff 1 (r). 

The left side of ( 1271) is an almost periodic function with respect to ou, in the Bohr's 
sense. Following the same arguments as those of the Theorem [1] we obtain 



3=1 \k^j 



Choosing ujj = vr/2rj and noting that, as the parameters {rj)^ =1 are not entire multiples 
of each other i.e. Ti/rj N: 



,71 Ti- 



cos(ujjTi) = cos( — — ) 7^ for i ^ j. 



2 Tj 



For each value ojj, we have 



] J cos 2 (-cj u.) / qj(s)ds = / qj(s)ds = 0. 

fc#j 'Jo Jo 

and finishes the proof of the Theorem [2J 
5. 5. Proof of Theorem 

In this subsection, we consider two potentials q = an d <? = ®f=iQ.ji satisfying 

the assumptions of the Theorem [5J Assuming that they give rise to the same Neumann 
reflection coefficients, R^{u) = Rtf(u), we have the characteristic equation (T23j) . 

Let us define the operator to be the operator —-£2+ qj(x) over [0, Tj] with the 
domain 

D(C^) = closure of C#(0, Tj) in H 2 (0, Tj ), 

where Cj^(0, Tj) denotes the space of infinitely differentiable functions / defined on [0, Tj] 
satisfying Dirichlet boundary condition at and Neumann boundary condition at Tj. 
By the perturbation theory for linear operators [11] . there exists e\ > small enough 
such that, if |M|z/>°(r), ||<?||l°°(T) < £i then the eigenvalues of £^ remain positive. 
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Considering ((A^) 2 )"^ (A| > 0) the sequence of eigenvalues of Cj^, (1231) implies for 
each j = 1, ■ • ■ , N, 

n^(0>M)^(0,Af) r%(x)<pb{x,)i)<&{x,>i)dx = 0, V< = 1,2,... (28) 

where we have applied the fact that ^^(0, A]) = 0. 

At this point, we will use the assumption Bl' on the lengthes Tj to obtain a Lemma 
on the non-overlapping of the eigenvalues for different branches: 

Lemma 3. Under the assumptions of the Theorem^ there exists e 2 > small enough 
such that, if \\q\\L°°(r), ||<?||L°°(r) < £2, then 

J] ^(0, \i)fa(0, K) ± 0, Vj = 1, • ■ ■ , N, V* = 1, 2, ■ • • . 

Mi 

Proof. In order to prove this Lemma, we only need to show that Xj is not an eigenvalue 
of £^y- nor for k ^ j. 

In this aim, we first show that there exists M > such that for integers ix,t2 > M, 
is different from Af and Af 2 the eigenvalues of and Assume, contrarily, that, 
for all M > 0, there exists z'i,Z2 > M, such that 

A^ = A^ or A£ = A*. (29) 

As the potentials g and g are if 1 , we have the following asymptotic formula's for the 
eigenvalues (see, for instance, Remark 1.1.1, page 7, j5]): 

A , = (2t - i)tt + /j^c^a + c /l 

2r fc 2ir fc \i 2 

- xk _ (2i-l)n | fft(# 1 O (- ) ■ 
2r fe 2ir fe U 2 



However, as by assumption B2, the integrals qj.(s)ds and qk{s)ds are zero, this 
implies: 

Therefore, the relation (129|) implies that, for all M > there exists > M such that 

Assuming, without loss of generality, that ii < i 2 and dividing the above inequality by 
(2ii — l)7r/27fc, we have 



(2ii - 1)tt (2i 2 - 1)tt 
2r 7 - 2r k 
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Therefore, we must have the existence of an infinite number of integer couples i 2 ) 
satisfying the above inequality. However, by Thue-Siegel-Roth Theorem [16], for any 
irrational algebraic number a, and for any 5 > 0, the inequality 



\a-p/q\ < l/\q 



2+6 

j 



has only a finite number of integer solutions p,q (q ^ 0). This, trivially, leads to a 
contradiction and therefore there exists M > such that for i±,i2 > M and j ^ k, 

Mi ± M 2 and Mi * Mr 

For the M first eigenvalues on each branch, we apply the perturbation theory for 
linear operators [TT]. Having j, the branch index, fixed, we will show that for e\ > 
small enough, if H^Hl^, H^IIl^ < e i then the quantities (\DfLi do note coincide with 
the quantities (Af)^ 1; where k ^ j. In fact, for the case of q = q = 0, this is, trivially, 
a consequence of the fact that the branch lengthes Tj and are 2-by-2 Q-linearly 
independent. Now, adding a small perturbation, q or q, in the generalized sense (see 
page 206 [IT]), by the continuity of a finite system of eigenvalues (see page 213 [H]), this 
claim remains valid. Thus, there exists a small enough e 2 = min J= i j ... > jv(e^) such that, 
if lklUoo> ll<?IUoo < e 2 then (A^)f£ 1 do not coincide with (Af)^ , for all j, k — 1, • • • , N 
verifying j ^ k. □ 

Applying Lemma [3] to the Equation (J28|) . for || g|| ioo , || g|| ioo < e := min(ei,e 2 ), we 
have: 



fy(xWjf(x, \l)<ffo(x, %)dx = 0, \/j = l,---,N, Vi = l,2, 
Using the integral representation with a Volterra kernel (1251) . we rewrite 

jT' ^cos 2 Xj (tj -x) + £ K^Tj-x^j-t) 2A ,-^ _ ^ \ ^ = 
^ ' ^-(jc) ^cos2A^'(T i - x) + y ' Tj ' - - cos2A|(r j - = 

So + / ^ 7 " J '~2 <,Tj '~ a?) ^^ d!C J cos2A K r i " *) dt 

where we have applied the assumption B2, qj = 0, for the passage from the second 
to the third line and the Fubini Theorem for the last passage. This implies that 



/ 3 A j (t) cos 2A|(r i - t)dt = Vi = 1, 2, 
Jo 



where 



o 
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By the assumption of Theorem [21 the potentials qj{t) are zero on [0, Tj/2] and 
consequently A^(t) = on [0, Tj/2]. Therefore, we have 

I 3 A 3 {t) cos2\ 3 (r j - t)dt = Vz = 1, 2, • • • . 

J Tj/2 

However, by the Proposition 1.8.6 of [5], the system {cos(2Aj(rj — .))}j=i,2,-- provides a 
Riesz basis over L 2 [rj/2, Tj]. Consequently, 

q.(t) + f '-^ J ' ] qM )<lr = t G [0, 

Jo * 

This homogenous Volterra integral equation has only the trivial solution qj = on [0, Tj]. 
This, trivially, implies that for e = min(ei,e 2 ), if |M|,L°°(r)> ||<?IU°°(r) < e then q = on 

r. 



5.^. Proof of Theorem^ 

We consider two potentials g = ®jLi9j an d <? = ®f =1 qj, satisfying the assumptions 
of the Theorem HI Assuming that they give rise to the same Neumann and Dirichlet 
reflection coefficients, Rj^(uj) = Rjv(uj) and Rx>(u) = Rd(u), we have the characteristic 
equations (1231) and (|2il). 

We define the operator C 3 V exactly as £jy (defined in the previous subsection) with 
Dirichlet boundary conditions at and at Tj . Still, by the perturbation theory for linear 
operators [TT] . there exists e\ > small enough such that, if ||qi|L°°(r), ||?|U°°(r) < e i 
then the eigenvalues of C 3 ^ and C 3 V are all positive. 

Considering ((X 3 ) 2 )Zi (K > 0) the sequence of eigenvalues of jC 3 j^, and ((/x^) 2 )^. 1 
(fi 3 > 0) the sequence of eigenvalues of C 3 V , the equations (T231 and ff2^1) imply: 

J] ^(0, >i)fa(0, A|) P <7,(*)^, Af W M {x, >()dx = 0, 

m .... Vj = l,---,iV,Vi = l,2,..., 

where we have applied the fact that (^(0, A]) = </4>(0> /4) = 0. Following the same 
arguments as those of the Lemma El there exists e 2 > small enough, such that 

k * j 3 = 1, • ■ • , N, % = 1, 2, • • • . 

n^>(o,/4)^,(o,/4)^o, 

Therefore, for each j — 1, • • • , iV we have , 

< &(0>A0> K) >L 2 (0,r J ) = < 9j(0>¥4>(-»^i) >i 2 (0,r,)= 0,Vi = 1, 2, • • • . (31) 
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At this point we define the sequence (v>i)j=i,-,N, t=o,i,2,- as follows 

^ (x) =(0,.-.,0, ^ ,0,..-,0), xGT 
j-th position 

4-i(*) = (o, • • • , 0, ^{x, \i)0ir(x, xi) - i 0, • ■ ■ , 0), x G r, i > 1 

V v ' 

i-th position 

4( x ) = (0, • • • , 0, <f? v (x, Hi)<Pv( x > A) ~ \> °> " " " > °)> x G T, i > 1 

V v ' 

j-th position 

As /J- 7 qj(s)ds = by Assumption B2, (jHTl) implies 

< ?(•), wj(-) >L 2( r)= 0, Vj = 1, ■ • • , N, Vi = 0, 1, 2, • • • . 

where 

N 

<V(-),W(-) > L 2 {T) = ^ < V j (-),W3(-) >L2(0,r J ) • 

i=i 

However, by Lemma [5] proved in the Appendix, the sequence (Wi)j=i,-,jv, i=o,i,2,- forms 
a Riesz basis over L 2 (T) and therefore g = 0, which finishes the proof of Theorem HJ 
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Appendix A. Riesz basis properties 

The goal of this appendix is to provide a Lemma [5] on the Riesz basis property for 
the sequences (t^i)j=i,—,N, i=o,i,2,~ defined in Subsection 15.41 Note that this is a direct 
consequence of a result already proved in [3]. However, we provide a proof for the sake 
of completeness. 

We first provide a classical result on Riesz sequences [SJ page 91]: 

Lemma 4. Let {fj}jL 1 be a Riesz basis for the Banach space B. Let {gj}'jL 1 be 
quadratically close to {fj}f = i, i.e. 

oo 

Then if the sequence {gj}jZ 1 is oo -linearly independent or if it is complete in B, then it 
is, also, a Riesz basis for B. 

Throughout the Appendix, we consider the sequence (ul) to be defined as in 
Subsection 15.41 We therefore have the following lemma: 
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Lemma 5. Assuming q and q in H 1 ^), the sequence (u J i )j = i t ... t N ! j=o,i,2,- provides a 
Riesz basis over the Hilbert space L 2 (T). 

Proof of Lemma{5[ Fixing the branch index j, we prove that {ul(x)}^ is a Riesz basis 
over L 2 (0,Tj). Here, for simplicity sakes, we have identified the vector u{(x) with its 
j-th component and therefore u\(x) is not a vector anymore but rather a function in 
L 2 (0, Tj). We also remove the branch index j. 

First, we prove that {vi(x) := f^-^u^x)}^ is a Riesz basis over L 2 (0,r). Note 

that 



T d t d 

Vi(0) = — -j-(<Pjv.2>)(x, \)<Pm,v(x, \) + — T-(<Ptf,v)(x, \)V>N,T>{x, Aj) 
m dx ' x=o m dx 

t d t d 

nr dx x=r m dx 



= 

x=0 



0. 



In fact, one can easily see that the functions are solutions of 

PX 

- v'l + 2{q{x) + q{x))vi + / N(x, s)vi(s)ds = 4A 2 ^, x E (0,r), (A.1) 



where 

iV(x, s) = g'(a;) + + (q(x) - q{x)) I (g(£) - g(0)d£. 



In particular, we obtain that the functions {vi(x)} c *i l are the eigenfunctions of (IA.ll) 
with boundary conditions Uj(0) = Vi(r) = 0. 

Then there exists a bi-orthonormal sequence {v*{x)}°^ l in L 2 [0,r], eigenfunctions 
of the adjoint operator 

-{v*)" + 2{q{x) + q{x))v* + / N{s,x)v*{s)ds = 4A 2 <, v*(0) = u*(r) = 0. 



Thus are cu-linearly independent. 

From the integral representation of fundamental solutions and from Lemma [U we have 

iirx ( 1 \ 

Vi(x) = sm( ) + O I - 1 as 2 ^ oo. 



So 

oo 



K(ar) -sin(^)|||2 (0jT) < oo 



E 

i=l 

i.e. {vi(x)}'?i 1 is quadratically close to {sin(— By the virtue of the Lemma HI 
this gives that {v^x)}^ is a Riesz basis over L 2 (0,r). 

Now let us show that the sequence of the functions {u^x)}^ is complete in L 2 (0, r). 
Indeed, for any / e L 2 (0,r), suppose that 

f(x)ui(x)dx = 0, i = 0, 1, 2, • • • . 
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In particular, considering it , we have 



f(x)dx = 0, 

Jo 

and therefore, 

/ f(x)(ui(x) + l/2)dx = 0, i = 0, 1,2, 
Integrating by part, we have 

\ X = T 

f(s)ds) (u t (x) + l/2) 



v i( x ) I / f( s )ds ) dx = 0, z = 1, 2, • • • . 



x=0 •'O 

and therefore using the boundary condition ■Uj(O) + 1/2 = 0, 

v i( x ) I / f( s )ds ) dx = 0, i = 1, 2, • • • 



The sequence {"^(x)}^ is complete, hence fj f(s)ds = for x G [0, r]. Therefore 
{■^(x)}^ is complete in L 2 (0,r). Since {r]iUi(x)}°l (where rj^s are appropriate 
normalizing constants) is quadratically close to {cos(^)}°i , it follows from Lemma 
HJthat {ui(x)}°Z is a Riesz basis for L 2 (0,r). □ 
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Abstract. The Frequency Domain Reflectometry (FDR) is studied as a powerful 
tool to detect hard or soft faults in star-shaped networks of nonuniform lossless trans- 
mission lines. Processing the FDR measurements leads to solve an inverse scattering 
problem for a Schrodinger operator on a star-shaped graph. Throughout this paper, 
we restrict ourselves to the case of minimal experimental setup corresponding to only 
one diagnostic port plug. First, by studying the asymptotic behavior of the reflection 
coefficient in the high-frequency limit, we prove the identifiability of the geometry of 
this star-shaped graph: the number of edges and their lengths. The proof being rather 
constructive, it provides a method to detect the hard faults in the network. Next, we 
study the potential identification problem by inverse scattering, noting that the poten- 
tials represent the inhomogeneities due to the soft faults in the network wirings. Here, 
the main result states that the measurement of two reflection coefficients, associated 
to two different sets of boundary conditions at the extremities of the tree, determines 
uniquely the potentials; it is a generalization of the theorem of the two boundary spec- 
tra on an interval [3] . 
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1. Introduction 

The rather extensive literature concerning the "inverse scattering problem" and the 
"inverse Sturm-Liouville problem" on graphs have mostly followed separate pathways 
except for a very few results [3, [6] . In the following paragraphs, we briefly recall the 
previous results on these subjects and at the end we will situate the result of this paper 
with respect to the others. Indeed, as it will be seen later, the inverse Sturm-Liouville 
problem considered in this paper raises from the necessity of finding a minimal setup 
for solving the inverse scattering problem. 

The paper [8] considers a star-shape graph consisting of N infinite branches and 
solves the inverse scattering problem assuming the measurement of N — 1 reflection 
coefficients. Next, in the paper [9], Harmer provides an extension of the previous 
result with general self-adjoint boundary conditions at the central node. This however 
necessitates the knowledge of N reflection coefficients. 

The paper [13] studies the relation between the scattering data and the topology of 
the graph. They show that the knowledge of the scattering matrix is not enough to 
determine uniquely the topological structure of a generic graph. In pQ, Avdonin and 
Kurasov consider again a star-shape graph with N finite branches. They, show that the 
knowledge of one diagonal element of the response operator allows one to reconstruct 
the potential on the edge corresponding to this element. 

As mentioned above, in parallel to the research on inverse scattering problems, 
another class of papers consider the inverse problem for Sturm-Liouville operators on 
compact graphs. These results can be seen as extensions of the classical result provided 
by Borg [3], on the recovering of the Sturm-Liouville operator from two spectra on a finite 
interval. The main progress in this field has been made by Yurko [IE] , [Hj and [20] . The 
paper [18] deals with the inverse spectral problem on a tree. The idea is to generalize the 
Borg's result in the following sense: for a tree with n boundary vertices, it is sufficient 
to know n spectra, corresponding to n different settings for boundary conditions at the 
extremities, to retrieve the potentials on the tree. In the recent paper [20], the same 
kind of result is proposed for a star-shape graph including a loop joined to the central 
node. Finally, [19J provides a generalization of [18] to higher order differential operators 
on a star-shape graph. 

In [15], the author proves that under some restrictive assumptions on the spectrum of 
a Sturm-Liouville operator on a star-shape graph with some fixed boundary conditions, 
the knowledge of this spectra can determine uniquely the Sturm-Liouille operator. 

In this paper, we consider the inverse scattering problem motivated by the 
application in fault-detection/diagnostic of star-shape LC transmission networks. We 
are interested in minimal experimental setting providing enough information on the 
network and the potentials on its branches. The graph consists of N finite branches 
joined at a central node and we add a infinite branch to this central node for the 
experimentation. We will see that the knowledge of the reflection coefficient is equivalent 
to the knowledge of the spectra for the Sturm-Liouville operator defined on the compact 
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part of the graph for various boundary conditions at the central node. By considering the 
associated inverse Sturm-Liouville problem, we will show that under some assumptions 
on the geometry of the metric graph, the knowledge of only two reflection coefficients, 
corresponding to two settings for the boundary conditions at the terminal nodes, is 
enough to determine the potentials (at least locally). 

In the next section, we explain the application under study and the associated 
experimental setup. We will show that the inverse scattering problem for the so- 
called Telegrapher's equation (lossless case) on the network is equivalent to an inverse 
scattering problem for a Schrodinger operator over the metric graph of the network. In 
Section [3l we will consider the direct scattering problem and we will characterize the 
reflection coefficient in terms of the fundamental solutions for Sturm-Liouville operators 
on branches. In Section HJ we will show that the knowledge of only one reflection 
coefficient is enough to identify the lengthes of the branches of the metric graph. This 
result will be useful to locate hard faults (open or short circuits). Finally, in Section [5] 
we consider the main problem of recovering the potentials from the knowledge of one or 
two reflection coefficients. This result will be useful to locate soft faults (local variations 
of the electrical characteristics). We will prove the equivalence of the inverse scattering 
problem with an associated inverse Sturm-Liouville problem on the compact part of 
the graph. This inverse problem is then treated applying the methods extending the 
classical result by Borg [3]. 

2. Frequency domain reflectometry 

The electric signal transmission through a wired network is, generally, modeled with 
the "Telegrapher's equation" and characterized by the parameters L, C, R, G (functions 
of the space position z along the transmission lines) representing the inductance L, 
capacitance C, resistance R and loss conductance G per unit of length. These parameters 
allow a rather complete and understandable description of the transmission lines and 
are sufficient to represent the lines in the frequency range used during reflectometry. In 
the sequel we will suppose that this model can be used for all frequencies. However, 
it appears to be impossible to retrieve all these parameters uniquely through the 
information provided by reflectometry experiments. Everywhere, through this paper, we 
will consider the simpler nonuniform lossless situation (R = G = 0). As we will see later, 
the reflectometry measurement is still not enough to retrieve the both parameters L and 
C but rather an aggregate of these two parameters, the local characteristic impedance 
Z c {z) := y/L(z)/C(z). 

Following [12] and [10], the presentation of the reflectometry experiment of 
this section, is based on a model derived from the "Telegrapher's equation" and 
parameterized by Z c . To cope with the network case, we have translated the Kirchhoff 
rules at the nodes of the network within this new modeling framework. Note that, in 
this paper, we restrict ourselves to the case of a simple star-shape network and therefore 
the only node of the graph where the Kirchhoff rules need to be adapted is the central 
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one. The faults, in which we are interested through this approach, are represented by 
the lengths of the branches (hard faults) and by the heterogeneities of Z c along the 
branches (soft faults). The considered reflectometry experiment model is based on a 
far-field method consisting in adding a uniform infinite wire joined to the network at its 
central node. In practice, connecting a matched charge to the extremity of a finite line, 
is sufficient to emulate the electrical propagation through an infinite line. 

The linearity of the transmission line model allows to replace any test by an 
equivalent test in harmonic regime. We can therefore start by stating the Telegrapher's 
equation in the harmonic regime, i.e. the tension and the current intensity are 
respectively of the form e~ tuJt V(u, z) and e~ tuJt I(u), z), where u is the time frequency 
and z the position. On each line, we have 

d 

—V(u), z) — iujL(z)I(uj, z) = 0, 

a (1) 

— I(u, z) - iuoC{z)V{uo, z) = 0, 

Through this paper, we assume that 

Al the distributed parameters C(z) and L(z) are twice continuously different iable on 
the transmission lines; 

A2 they are strictly positive, C(z) > 0,L(z) > 0; 



A3 the characteristic impedance Z c (z) := yL{z)/C{z) is continuous at the central 
node of the star-shape network; 

A4 the transmission lines are uniform in a neighborhood of the extremities of the 
branches. 

The Liouville transformation. Note that the reflectometry experiment leads to 
observing the tensions and currents along the time at some position: only the travelling 
times (and amplitudes) of waves are accessible by such experiments. A fault can only 
be localized in terms of the traveling time of the reflected test wave starting from the 
test point. This leads to a particular change of variables, the Liouville transformation, 
allowing to work with the traveling time rather than spatial coordinates. Let us recall 
this transformation: 

x(z) = [ y / L(z)C(z)ds 
Jo 

which corresponds to the wave traveling time from the position to the position z. 
Remark that after this transformation, u is also the wave number on each branch. 

The inverse transformation being well defined, we will write C{x) = C(z(x)), 
L(x) = L(z(x)), V(u,x) = V(u),z(x)) and I(u,x) = I(u,z(x)). 

The Telegrapher's equation (TjQ) becomes 

d 

— — V(u, x) — iuZ c (x)I(u>, x) = 0, 

d , (2) 

— x) — wjZ c (x) 1 V(lu, x) = 0. 

(y JU 
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The wave decomposition and equivalent forms of the Telegrapher's equation. Define 

1 v ( ^_ x dZ c {x) 

dx 

following decomposition of V: 



V±{uj,x) = \{V{w,x) ± Z c (x)I(u,x)) and Q{x) = -Z c {x)- 1 £ ; . We have the 



V(u, x) = V + (uj, x) + VL.(w, x), 

d _ (3) 

— V±(u, x) =F vu)V±(u, x) = ±Q(x)(V + (uj, x) - V-(u, x)) 

In particular, in an interval where a branch is uniform, Q(x) =0, and the solution 
is the sum of waves of opposite directions. For any x and x in this interval: 

V(u,x) = V + (to,x)e tuj{x - s) + V„((j,x)e- tuj{x - s) . 

Define now y(u>,x) = ^Y c (x)V(cu, x) and q(x) = y/Zjx)-^ y/YJjc) , with Y c {x) = 
Z c (x)~ l . The Telegrapher's equation (CQ) becomes a Schrodinger equation: 

d 2 y 

- -j^(uj,x) +q(x)y(uj,x) = u 2 y(u,x). (4) 

It can be seen from (jHJ) and (jlj), that, the knowledge of the potential Q(x) or of q(x) 
and of the boundary conditions on V and I is sufficient to compute the solution of (PQ) on 
the network. In our lossless situation, we have chosen q as the parameter to be identified 
through the reflectometry experiment. A variant of ([3]), the Zakharov-Shabat equations, 

would be the good choice in the more general lossy case. Remark that — Q 2 + q = 0. 

dx 

The reflection coefficient. Taking e tu ( x ~ 1i \ oriented toward the increasing x, as 
a reference forward wave, the reflection coefficient is the following ratio of backward 
over forward wave amplitudes: R(u,x) = v~(ux)e-*» x = e 2 ^ x ' ^~|^ . In particular, 
R(u,x) = R(u,x): defined in this way, the reflection coefficient is constant in intervals 
where Z c is constant. For an arbitrary x, denoting by Z(u,x) = V(u,x)/I(u,x) the 
(possibly infinite) apparent impedance at x, we still define the reflection coefficient as 
being 

RM=e^ T X) - Z 7 fy (5) 
Z(uj, x) + Z c (x) 

With this definition, in general, V-(u,x) = e~ 2tuJX R(u,x)V + (uj,x), and, using Q, it is 
easy to check that R is solution of the following Riccati equation: 

9B &> X ) - e - 2wx Q(x)R(LU, x) 2 + Q{x)e 2wx = 0. (6) 
ox 

Finally, note that, we will only consider "positive real" terminal impedances Z(u,t), in 
the sense that Z(— u, r) = Z(u,r) and %tZ(u), r) > 0. This together with the fact that 
Z c [x) G M. is positive implies that \R(uj,x)\ < 1. Furthermore, experiments with uj > 
are sufficient. 

The network under test. Throughout this paper, F represent the compact star- 
shape network consisting of the branches {ej)^ =l joining at the central node and r + is 
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the extended graph where the test branch e is also added to the graph. We have N + 1 
equations of the form 

d 2 y 

- + Qj(x)yj = u?yj x G (0, Tj), (7) 

where Tj is the wave traveling time associated to the branch number j (to = oo 
as the added branch eo is assumed to be an infinite line). In particular note that, 
as the infinite branch eo is assumed to be a uniform transmission line, we have 
qo(x) =0, x E (0, oo). 

The boundary condition for the reflectometer. Consider now that a generator with 
a matched internal impedance Z c and an electromotive force 2V g e tuJX is connected in x of 
an interval X where Z c is constant. We have V(lu,x) + Z c (x)I(lu,x) = 2V g e tulx which can 
be the boundary condition if the branch terminates at x. We have V + (u,x) = V g e lu)X 
and V(u,x) = V g {e tuJX + R{u,x)e~ tuJX ) for all x e 1 . As, R(u,x) = R(u,x), the 
reflection coefficient can then be determined from the measurement of the tension 
V(uj,x) anywhere in X. In the sequel, we will use a test branch eo with a constant 
Z c connecting a matched generator to the central node . The measured reflection 
coefficient on this branch, will be simply written R(u). As we will work with y, we 
choose V g = a/ Z c (x), so that y(uj,x) = e luJX + R(uo)e~ tulx . Finally, it will be convenient 
to take the same positive orientation on all the branches, from the central node at x = 
toward the increasing x. Our reference forward wave on e is then in the direction of 
the decreasing x, so that, changing x into —x, and supposing e of infinite length, the 
boundary condition for the reflectometer is: 

y(u, x) ~ e~ luJX + R(u)e wx as x -> oo in e . (8) 

The two sets of boundary conditions at the network extremities. In order to recover 
the potential of the star-shape network, we will need to consider two experimental 
settings, with open circuit or short circuit at the extremities of the branches. This will 
lead to a problem similar to solving an inverse spectral problem for the Sturm-Liouville 
operator when two spectra are known. 

The first setting corresponds to open circuit configuration at the extremities of the 
finite branches ((e-,)^). This, together with the Assumption A4 on the local uniformity 
of the lines around r/s, leads to boundary conditions of the form Ij(cu,Tj) = 0, or 
equivalently, we obtain the setting called, the Neumann configuration: 

2/;.(r,) = j = l,...,N. (9) 

The second setting corresponds to the short circuit configuration at the extremities 
of the finite branches ((e.,)^). This leads to boundary conditions of the form 
Vj(u,Tj) = 0, or equivalently, we obtain the setting called, the Dirichlet configuration: 

y j (T j ) = J = 1,---,AT (10) 

Remark 1. In some of the applications that we have in mind, the reflectometry 
experiment has to take place without perturbing significantly the normal utilization 
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of the transmission network, so that using open or short circuits conditions would be 
impossible. There is a way to circumvent this problem by computing the results of the 
open or short circuit experiments from results of two less invasive experiments. The 
idea is to use nonlinear superposition properties of solutions of Riccati equations as in 
[17] . in order to get a closed-form representation of the reflection coefficient, solution of 
(jHJ), as a function of a general load impedance (value of Z at the extremity of a branch) 
and of two particular solutions corresponding to two load impedances more compatible 
with the network utilization. This will be presented in a forthcoming paper. 

The boundary condition at the central node. It writes 

yi(u,0) = Uj(u,0) =: y(u) i, j = 0, • • • , N, 

2^2/i(w,0)-2/ (w,0) = -- z0 

j=l c 

where y'j (u,0) and (Z^')'(O) denote the spatial derivatives at the point x = and Z 3 C is 
the characteristic impedance of the branch number j. Note, in particular, that we have 
applied the continuity of Z^s at the central node (Assumption A3): Z 3 c {jS) = Z®, Wj. 

Formulation of the model. In conclusion, in order to study the LC-transmission 
line equations on the graph T + , we can study the Schrodinger operators 

d 2 

D(C^ V ) = closure of C% v in H 2 (T + ), (12) 

where C^(T + ) (resp. C^(T + )) denotes the space of infinitely different iable functions 
/ = ®^=ofj defined on Y + satisfying the boundary conditions 

./•«') ./••(()) j,/ = 0,---,A, 



E ft® - /o(0) = Hf (0), H = -- V zQ L (13) 

j=l c 

/j(ri) = (/ J (r J )=0for^(r + )), j = !,•••, A. 



3. Direct scattering problem 

The operators (C^ v , D(Cjj- v )) are essentially self-adjoint. To prove this fact we observe 
first that these operators are a compact perturbation of the operators ®" =0 ( — j with 
the same boundary conditions. Now, we apply a general result by Carlson [4J on the self- 
adjointness of differential operators on graphs. Indeed, following the Theorem 3.4 of pE], 
we only need to show that at a node connecting m edges, we have m linearly independent 
linear boundary conditions. At the terminal nodes of {ej}^ =l this is trivially the 
there is one branch and one boundary condition (Dirichlet or Neumann). At the central 
node it is not hard to verify that (llip define iV + 1 linearly independent boundary 
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conditions as well. This implies that the operators (C^- v , D(Cj^ v )) are essentially self- 
adjoint and therefore that they admit a unique self-adjoint extension on L 2 (T + ). 
We are interested in the scattering solution where a signal of frequency uj is applied at 
the infinite extremity of the infinite branch. In such a case, we will be seeking a solution 
satisfying the asymptotic behavior 

y (x, uj) ~ e~ luJX + R(uj)e l " x , for x -> oo. 

Proposition 1. Under the assumptions Al through A4, there exists a unique solution 

continuous with respect to uj, of the scattering problem, satisfying 

• ~~£zyj( x i u ) + 1j( x )yj( x ^ w ) = uPVjfa for j = 0, ■ ■ • , N; 

• (yj(x,uj))j =0 satisfy the boundary conditions (Tl\) and (G|) or ( T70j) ; 

• For each uj G R, there exist R(uj) such that 

y {x, uj) ~ e~ luJX + R{uj)e iujx , x -> oo. (14) 

VKe wz// denote the reflection coefficient R(uj) defined by (fl^[ ) m i/ie Neumann (resp. 
Dirichlet) case by Rj\f(uj) (resp. Rd{uj)). This coefficient appears to be unique. 

Proof. This proof gives us a concrete method for obtaining scattering solutions. Indeed, 
we will propose a solution and we will show that it is the unique one. 
In this aim, we need to use Dirichlet /Neumann fundamental solutions of a Sturm- 
Liouville boundary problem. 

Definition 1. Consider the potentials q 3 as before and extend them by on (— oo,0) 
so that they are defined on the intervals (—00,7^]. The Dirichlet (resp. Neumann) 
fundamental solution (p J v (x,uj) (resp. <p J j\f(x, uj)), is a solution of the equation, 

d 2 



rVW x ' <*>) + q^x^Jx, uj) = uj f^Jx, uj), x G 



"OO, T; 



^ v {Tj,Uj) = 0, (^)' v (Tj,Uj) = 1, 

vm^j,^) = i, <y )V( r i> u ) = °- 

Consider, now, the function 



where 



$! Vj jv(x,oj) = ® j=o y J vM (x,u), 
^° vM (x, uj) = e~ tuJX + R vM (uj)e wx , x G [0, 00), 
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Here the coefficients Rv,jv and o? v ^ are given by the boundary conditions (1111) at 
the central node: 

l + R Vf M-(u) = c4) >Af {uj)(fi)^(0,Lo), j = l,---,N, (15) 

N 

Y J ^ vM {u){i ( P v ^)'{Q,oj) + - ife^M) = H(1 + Rdst(u)). (16) 

3=1 

One, easily, sees that this ^v,Af is in D(£j^ v ), the domain of the operator, and satisfies 

the conditions of the proposition as soon as the coefficients Rv,M and {pSpjJ) are 
continuous with respect to u>. This, trivially, provides the existence of a scattering 
solution. Here, we show that ^t>,jv is actually the unique one. 

Assume that there exists another Y-djj = ®jLo^r>,Ar( x i u ) solution of the scattering 
problem. By the compact injection of H 2 in C 1 , we now that Yp^x, u) and ^/ VJ ^(x, u) 
are C 1 functions of x over [0, 7}]. Therefore, the Wronskian 

is well-defined. Moreover, as the potentials qj(x) are continuous functions over [0, Tj] and 
as Ypfj-(.,uj) and ^ 3 t> k[{-, are solutions of the associated Sturm-Liouville equation, 
they are in fact C 2 over [0, Tj]. Thus, the derivative of the Wronskian is also well defined 
over [0, Tj]. Through a simple computation and by noting that Yp u>) and ^ 3 T) u-{.,uj) 
are solutions of the same Sturm-Liouville equations, one has 

-^W(Yi M (x, u), %^{x, u)) = 0, x G [0, Tj ], 

and so the Wronskian remains constant over the interval [0, Tj]. 

For the finite branches {ej)^ =1 , applying the (Dirichlet or Neumann) boundary 
conditions at the terminal nodes, we easily have 

W{Yl^{ T j,uj),¥ ViM {Tj,u)) = Q, 

and therefore the Wronskian is identically over the whole branch. This implies that, 
Y*(.,uj) and % ^(., 10) are co-linear: 

Y vM (x, u) = (3 J v ^{uj)(f vM (x, u), x e [0, Tj], j = 1, • • • , N. 

Over the branch eo, as Ypj^(.,uj) satisfies a homogenous Sturm-Liouville equation 
{Qo = 0); it necessarily admits the following form 

What remains to be shown is that one necessarily has Rv,Af{u) = Rv,Af{^) and similarly 

Indeed, the equations (fT5j) and ffl6l) provide N + 1 linear relations for the N + 1 
unknown coefficients Rv,m and Trivially, as soon as, the coefficients 
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{ ( Pv at(^> u; ))j-i are non " zero ! these linear relations are independent and there exists 
a unique solution for the unknowns R-pjV and («^, A ^)^_ 1 . However, the zeros of each 
one of the coefficients {tfP v ^(0, correspond to isolated values of uj (square-root of 

the eigenvalues of the operator — +qj(x) with Dirichlet boundary condition at x = 
and Dirichlet or Neumann boundary condition at x = tj). Therefore, the coefficients 
Rv,m and (o4?,Af)^-i are well-defined except for a set of isolated values of uj. However, 
as these coefficients need to be continuous with respect to u), they will, also, be defined 
uniquely over these singular points. 

Furthermore, dividing (|T5|) by {1 + R v ^{uj)) and inserting (1T5|) . we find the explicit 
formula 

1 - RvMu) _H 1 A (^)'(O, 
Finally, inserting the value of Rd^{uS) found in (|T7|) into (fT5|) . we find 



□ 



4. Inverse scattering for geometry identification 

As a first inverse problem, we consider the inversion of the geometry of the network. In 
fact, we will prove the well-posedness of the inverse problem of finding the number of 
branches N and the lengthes (Tj)^ =1 of a star-shape graph through only one reflection 
coefficient Rj^(uj) (the case of Dirichlet reflection coefficient can be treated similarly). 
Furthermore, as we will see through the proof of the Theorem [TJ the method is 
rather constructive and one can think of an algorithm to identify the lengthes, at 
least approximately. The proof is based on an asymptotic analysis in high-frequency 
regime of the reflection coefficient and some classical results from the theory of almost 
periodic functions (in Bohr sense). Before, announcing the main Theorem, we need 
a few lemmas. A first lemma precises the high frequency behavior of the Neumann 
fundamental solutions (yr^-)^: 

Lemma 1. Consider a potential q in C°((— oo, r]) and take the Neumann fundamental 
solution, (pj\f(x,u), defined as in Definition^ We have 

<pj\f(0,uj) = cos(cjt) + £>(-), as uj — > oo, 

w 

(<pj\f)'(0, uj) = u sin(a;r) + O (1), as uj — > oo, 
where (<f_\f)'(0, uj) denotes the spatial derivative 4-cpj\f(x,uj) at x = 0. 
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Proof. We start by writing <fjy(x, u) in its integral representation. Indeed, the 
fundamental solution (pj^(x, uj) is given by [H]: 

(prt(x, uj) = cos(lu(t — x)) + / Gj^{t — x, r — s) cos(u;(t — s))ds, 

J x 

where G_\f(x,y) is a real function with the same regularity as jj q(s)ds: here, it is C 1 
with respect to both coordinates. We note that, as Gj\f is in C 1 , by integrating by parts, 
one has: 

r 1 s=t 

/ Gj^(t, t — s) cos(u;(r — s))ds = GV(t, t — s) sin(u;(r — s)) 

Jo u 3=0 

+- f T T-( G ^(^ T - s )) sin (^ - s))ds = o(-\ 
uj J ds \ujJ 

Therefore ipj\f(Q,uj) = cos(ujt) + O (1/co) as u> — > oo and we have the first relation. 
For the spatial derivative (<pf\f) r (x,uj) at the point x = 0, we have: 

((pjf)'(0, uj) = uj sin(c<j(r))— G(r, r) cos(u;r)+ / — [Gj^ij— x, r— s)) cos(u;(r— s))ds. 



dx x=o 



The kernel y) being C 1 , we have 



r d 

— G(t, t) cos{ujt) + / — (Gjs[{t — x, t — s)) cos(u(t — s))ds — O (1) as UJ —>■ oo. 
Jo x=o 

Thus, (</Pj^)'(0, a;) = u;sin(a;(r — x)) + O (1) as uj — > oo and the second relation 
follows. □ 

As we see by Lemma UJ in the high-frequency regime, the fundamental solutions 
become asymptotically independent of the potential q. The next lemma, provides an 
explicit method to identify the number N and the lengthes {tj)^ =1 of the branches for 
the homogenous case, where we know that the potentials {qj)f = i are all zero. 

Lemma 2. Consider a star-shape network T composed of rij branches of length Tj 
(j = 1, • ■ ■ ,m) all joining at a central node so that the whole number of branches N 
is given by Y^j=i n r Assume the potential q on the network to be (q = 0). Then 
the knowledge of the Neumann reflection coefficient Rj^{uj) determines uniquely the 
parameters {nj)^L 1 and (t^)^. 

Proof. We need to apply the explicit computation of the reflection coefficient provided 
by (fTTl) . The fundamental solutions are given, simply, by ipj^(x,uj) = cos(u;(r 7 - — x)). 
Therefore: 

1 — R\f{uj) 1 rT 1 v-^ ujsvairjuj) 

-H — y nj — 



1 + Rj\[(uj) iuj iu cos(rju;) 
The knowledge of Rm{uj) determines uniquely the signal: 



3=1 
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Assuming, without loss of generality, that the lengthes Tj are ordered increasingly 
7"i < • • • < r m , the first pole of the function f(u>) coincides with 7r/2r m and therefore 
determines r m . Furthermore, 

n m = lim cos(W m )/(u;), 

and therefore one can also determine n m . Now, considering the new signal g{oj) = 
f(u>) — n m tan(co>r m ), one removes the branches of length r m and exactly in the same 
manner, one can determine r m -\ and n m -\. The proof of the lemma follows then by a 
simple induction. □ 

We are now ready to announce the main theorem of this section: 

Theorem 1. Consider a star-shape network T composed of nj branches of length Tj 
(j = l,---,m) all joining at a central node so that the whole number of branches 
N is given by ^2J l = i n j- Assume the potential q on the network to be, simply, C°. 
Then the knowledge of the Neumann reflection coefficient Rj\f(uj) determines uniquely 
the parameters (nj)™^ and (Tj) 1 JL 1 . 

Proof. Assume that, there exists two graph settings (r J -,gj)^ 1 and {fj,qj)^ =l (the 
lengthes Tj are not necessarily different) giving rise to the same Neumann reflection 
coefficients: Rj^(cu) = Rj^{oj). By the explicit formula f|T71) . we have 

V (^)'(0, W ) _ if(^)'(o )W ) 



(,i ^—^ ,J Yn ,.,\ 1,1 j—~< 



This is equivalent to: 

JV / N 

n^(o,w) ^(^y(o^)n^(o^) 

3=1 \k=l tyk 

N I N \ 

n ^(o, U ) ^(^)'(o, u) n ^(o, u) = o. (is) 

3=1 \k=l tyk J 

Defining the function: 

N / N \ N I N \ 

F(oj) := Y\ cos(wf fc ) ^sin(cjr fe ) JJcos(urrz) - JJcos(an-j-) ^ sin(wf fc ) cos(a;f z ) , 

3=1 \k=l l^k ) j=l \k=l lj=k ) 

the asymptotic formulas of Lemma [1] and the ffl8l) imply 

F(uj) = 0{1/uj) ascj^oo. 

However, the function F(u) is a trigonometric polynomial and almost periodic in the 
Bohr's sense [2]. The function F 2 (u) is, also, almost periodic and furthermore, we have 

M(F 2 ) := lim ^ / F 2 {u)du = lim ^ ( I F 2 {u)du + [ F 2 {u)du 

Sl^ooiZ Jq C-^oo S I \Jq Ji 

1 „ ] 



< lim - d + C 2 / — du; ) = 0. 
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This, trivially, implies that F = (one only needs to apply the Parseval's Theorem 
to the generalized fourier series of the function F). However, the relation F(u) = is 
equivalent to 



N N 



i=i i=i 

and therefore, by Lemma [2], the two settings are equivalent and the theorem follows. □ 



5. Inverse Scattering for potential identification 

A second inverse problem, related to the detection of soft faults in the network, can be 
formulated as the identification of the potentials on the branches. Here, we consider 
the case of homogenous perfect wires. We will show that the measurement of the two 
reflection coefficients Rv{oj) and Rjyiu), corresponding, respectively, to a short circuit 
and an open circuit experiment, is enough to identify uniquely the small changes in 
the potential remaining in a certain regularity class. Indeed, we will prove that the 
inverse problem of finding the H 1 ^) potentials in an e L°°(r)-neighborhood of the zero 
potential (homogenous case), is well-posed for e small enough. In this aim, we will need 
an additive assumption on the electrical lengths (rj)j =1 , in order to remove symmetries 
leading to degeneracy problems. However, as it will be discussed later, it seems that 
this assumption can be relaxed. 

In a first result (Theorem [2]), we prove that under some natural assumption on 
the electrical lengths, the knowledge of only one reflection coefficient (here Rj^(u)) is 
sufficient to identify uniquely the values: 

Qj = / Qj{s)ds Vj = 1,---,N. 
Jo 

If anyone of these quantities appear to be different from zero, we know that there has 
been a change of parameters in the corresponding branch. By performing classical 
inverse scattering techniques over this branch [T2l [TO] we can identify its soft faults. 
The question is therefore to identify the soft faults in the network which does not 
change the quantities (qj)jL v This is treated in the two Theorems H] and [3j For both 
these theorems, we will need some more restrictive assumptions on the electrical lengths 
{ T j)jLi- I n Theorem El we will show that the knowledge of only one reflection coefficient 
(here Rjy(u)) is enough to identify the potential when we know that it has not changed 
on the first half of the branches. The Theorem H] provides a well-posedness result for 
the inversion of the potential over the whole graph but necessitates the knowledge of 
both reflection coefficients Rx>{oj) and R^{uj). 

In the sequel, we note that the potential over the infinite branch eo is always given 
and is 0. Indeed, this homogenous line is added only for the reflectometry experiment. 
The following theorem provides a global inversion result concerning the quantities q~j. 

Theorem 2. Consider a star-shaped graph Y and assume that 
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Bl The electrical lengths {rj}^ =1 are not entire multiples of each other, 

j- ■ 

— £ N for i^ j. 

n 

If there exist two potentials q = ®^ =1 qj and q = ®^ =1 qj in H 1 ^) giving rise to the same 
reflection coefficient, Rj^ioj) = Rj\t{uj), one necessarily has: 

qj(s)ds= qj(s)ds j = l,---,N. 
Jo 

Remark 2. Note that the result of the Theorem is also valid if we have Rt>{oj) = i?x>(a>). 

This theorem allows us to identify the situations where the soft fault in the network 
causes a change of the quantities q~j. As explained above, the inverse problem can 
then be considered on separate branches and solved through classical inverse scattering 
techniques. In the two following theorems, we assume that the soft faults in the network 
leave the quantities q~j unchanged and as for the perfect situation we are dealing with 
homogenous lines we will assume that 

B2 qj = qj (s)ds = for j = 1, • • • , N. 

Theorem 3. Consider a star-shaped graph F and assume that 

Bl' For any j,j' G {1, • • • , iV} such that j ^ j' , Tjjry is an algebraic irrational number. 

Then, there exists e > small enough such that, if the potentials q and q belonging to 
and satisfying B2, ||<2i|z,°°(r), ||<?IU°°(r) < e and qj(s) = qj(s) for s e [0, Tj/2], 
j = 1, • • • , N, give rise to the same Neumann reflection coefficients, Rj^(cu) = Rj^{uj), 
then q = q. 

Remark 3. Once again, the result remains valid if we replace the Neumann reflection 
coefficient by the Dirichlet one. 

Remark 4. The assumption Bl' seems very restrictive. However, as it will be seen 
through the proof, the only thing we need is that for any of fractions tj/tj/, there exists 
at most a finite number of co-prime factors (p, q) E N x N, such that the Diophantine 
approximation 

IL _ P 
Tj> q 

holds true. However, this is a classical result of the Borel-Cantelli Lemma that for almost 
all (with respect to Lebesgue measure) positive real a's this Diophantine approximation 
has finite number of solutions. 

Theorem 4. Consider a star-shaped graph F satisfying Bl'. There exists e > 
small enough such that, if the potentials q and q belonging to H 1 (F) and satisfying 
B2 and ||qi|L°°(r), ||?IU°°(r) < e, give rise to the same Neumann and Dirichlet reflection 
coefficients, 

Rj\f(uj) = RjKf(ou) and Rv{w) = Rt>(u), 

then q = q. 



< 
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Remark 5. The Theorem H] is a natural generalization to the case of a graph of the 
Theorem of two boundary spectra on an interval [3] (see, for instance, Theorem 1.4.4, 
Page 24, [5]). In [3j, Borg proved that the knowledge of two spectral data corresponding 
to two boundary conditions, determine uniquely the potential on an interval. Here 
the reflection coefficients Rj\f(u) and Rv(uj) play the role of this spectral data (see the 
Subsection 15. ip . 



5.1. Inverse Sturm- Liouville problem 

Throughout this subsection, we will consider a general star-shaped graph V (of N finite 
branches) and a potential q = ®f =1 qj belonging to H l (T). Furthermore we assume for 
the potential q that, the norm ||g||L°°(r) is sufficiently small. 

The main objective of this subsection is to show that the knowledge of the reflection 
coefficient R/^(uj) for C^- (resp. Rv(u) for £j) is equivalent to the knowledge of 
different positive spectra of Sturm-Liouville operators defined on T with Neumann (resp. 
Dirichlet) boundary conditions at terminal nodes and for various boundary conditions 
at the central node. In fact, defining the function 



[1 + Rjw(w)) ' 
where H is given by ( 1131) . we have the following result. 

Proposition 2. Fix w G I and define the Schrddinger operators Ctfp(uj) on the 
compact graph V as follows: 

d 2 

D{C M>v {u)) = closure of C^ V (T) m H 2 (T), 

where C^?^(r) (resp. C^ V (T)) denotes the space of infinitely differentiate functions 
f — ®f=ifj defined on T satisfying the boundary conditions 

/;(0) //(0) ■■/ j,j' = l,---,N, 

N 

j'=i 

/ifo) =0 (/ J -(r i ) = 0/ O rC^,(T)), J = l,---,N. 

Then we are able to characterize the positive spectrum of Cn^u) as a level set of the 
function h^^{oj): 

Proof. We prove the proposition for the case of Neumann boundary conditions. The 
Dirichlet case can be treated exactly in the same manner. We start by proving the 
inclusion 

o- + (C x (u)) C {e | f G R, MO = M«>)} • 
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Let £ 2 G a + (C^(u)), then there exists \I/ eigenfunction of the operator Cj^{ui) associated 
to £ 2 . In particular, it satisfies 

N 

2*;.(o) = m<*>)*, 

where \I/ is the common value of $ at the central node. 

Now we extend ^ to the extended graph r + , such that is a scattering solution 
for Cjj (see the Proposition [1]). In particular, the function must satisfy, at the 
central node, 

*t(0) = *+(0), j = l,---,N, 

N 

^(^)'(0)-(^ + )'(0)=^+(0). 

i=i 

Noting that \l/ is an eigenfunction of (C_\f(uj), D(£n{uj)), we have 

TV 

M"W(0) - W)'(O) = - W) = (0). (19) 

Now, noting that over the infinite branch admits the following form 

= R M {i)e lix + e~* x x G [0, oo), 

the relation ( fl9l) yields to 

/wHOMO + 1) - <(MO - 1) = #(M0 + 1), 

or equivalently 

This proves the first inclusion. Now, we prove that 

<t + (£a/-M) 3{e 2 |eeR, MO = M")} • 

Let ( G 1 be such that TiaKO = We consider a scattering solution of 

the extended operator (defined by ([ED ) associated to the frequency £ 2 . We, then, 
prove that the restriction of to the compact graph V is an eigenfunction of Cj^{uS) 
associated to the eigenvalue £ 2 . This trivially implies that £ 2 G a + (Cj K f(u)). 

In this aim, we only need to show that this restriction of \& + to T is in the domain 
D{C^{uj)). Indeed, this is equivalent to proving that the boundary condition: 

JV 

5>+)'(0) = /^)*o + (0), (20) 
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5>+)'(0) = H*+(0) + (*+)'(0) = (if + (0). 



Furthermore, 



tt+(0) = Rjsit) + 1 and (*+)'(0) = t^AO ~ 1) 



and so 



JV 



*+(o) = MO«tf(o) = /iH^(o). 



□ 



Rm(0 + 1 

This proves (1201) and finishes the proof of the proposition. 

We have shown that, the knowledge of the reflection coefficient Rtfp(uj) for T + is 
equivalent to the knowledge of the positive part of the spectra for £/v",d f° r ah values h 
at the central node. 

The following proposition provides the characteristic equation permitting to identify 
the eigenvalues of the operator Cj\f jT> (uj): 

Proposition 3. The real A 2 > is an eigenvalue of the operator Ctf t v(w) if and only if 
^at,»(A) = fy/vr,x>M$A/VD(A), 

where 



i=i 



i=0 



(21) 



(p^f£>(x, A) feeing t/ie fundamental solutions on different branches. 



Proof. We give the proof for the Neumann boundary conditions, noting that the 
Dirichlet case can be treated, exactly, in the same manner. Assume A 2 to be a positive 
eigenvalue of C^{u). The associated eigenfunction, y\{x) = has necessarily 

the following form: 

where a/s are real constants and the vector (ax, ■ ■ ■ , ajv) is different from zero. The 
function y x , being in the domain D(Cj\f(u)), it should satisfy the associated boundary 
condition at the central node. This implies that the vector • • • , a^) is in the kernel 
of the matrix: 



M 



Mr(0,A) 





-<^(0,A) 
<^(0, A) -<^(0, A) 
<^(0, A) 



\-Mw)y>V(0, A) + ^(0, A) ^2,(0, A) ViV(0, A) 



\ 



W0, A) / 
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where ^-(0, A) denotes -^<pA x ' A) \ x =o- This means that the determinant det(M) is 
necessarily 0. Developing this determinant we find: 

^(A) = /wM«MA). 

□ 

Corollary 1. Consider two potentials q = <%>f=iqj and q = (8>jLi9j and denote by Cj^- 
and the associated Neumann Schrodinger operators defined on the extended graph 
T + . Assuming that the of the reflection coefficients Rj^{uS) and Rm{oj) are equivalent 
Rj\f(u>) = Rtf(u), we have 

$/tf{w)Vtf{u)) = $jv(w)*jv(w), Wu e R, (22) 

where \l/_/v-, ^JV an d are defined through (j2*T|) for the potentials q and q. 

Proof. By Proposition [21 uj 2 is an eigenvalue of the operator Cm{uj) and £jv(a>). 
Applying the Proposition [31 this means that 

^f(u) = hx(ui)<&tf{u) and ^!{uj) = hj^(u)^j^{Lo). 

As Rjj(uj) = Rj\f(uj), we have hj^(ou) = hj^iu) and thus the above equation yields 
to (El. □ 



The above corollary is also valid when we replace the Neumann by Dirichlet 
boundary conditions. Finally, this corollary yields to the following proposition on the 
difference between the two potentials q and q. 

Proposition 4. Consider two potentials q = ®jLi?j and q = (8>jLi?j and denote by Cjj 
and Ctf, the associated Neumann Schrodinger operators defined on the extended graph 
T + . Assuming that the of the reflection coefficients R_\f(u>) and Rj^{uj) are equivalent 
Rtf{yj) = Rj^iyj), we have 

N m 

X)II^ ' W ^W ' W ) / ^A x WA x ^WA x ^) dx ^^ Vwgk, (23) 

3=1 k+j J ° 

where qj = q~j — qj . 

Proof. For j — 1, • • • , N, we have: 

q(x)<f^ r (x,u)(p^ r (x,u)dx- / qj(x)(ffa(x,w)(fo(x,w)dx = 



JO 

= vA x ' u )^A x > u )\ x x=o - ^A( x ^)A( x ^)\ x x=o 

= ^(0, utfAo, u) ~ ^(0, w)^.(0, w). 
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Here the second line has been obtained from the first one, replacing qj(x)(p 3 J y(x,u) by 
-^cfp^x, to) + to 2 Lpj^(x, to) and integrating by parts. Using fl22|) and the above equation, 
we have: 

N p 

3=1 k+j J ° 

N 

II ¥&(°> ^(°> w ) (^(0, w)<^(0, u) - ^(0, ^)^(0, a;)) = 
j=i k^j 

^^(to)®^) - $ M (to)^^(to) = 0. 

□ 

Before finishing this subsection, note that, once more, the above proposition is also 
valid for the case of Dirichlet boundary conditions and Rt>(uj) = Rx>{to) implies: 

N m 

22Yl ( Pv( Q ^)'fU ^) / qj(x)<p } v (x 1 u)0' 13 (x 1 u)dx = O, \/to e R. (24) 
We are now ready to prove the Theorems [2} [3] and HI 
5.2. Proof of Theorem^ 

We prove the Theorem [2] applying the characteristic equation ( l23|) and the integral 
representation of the functions iffo v (x , uj) . Again, for simplicity sakes, we give the 
proof only for the case of Neumann boundary conditions, noting that the Dirichlet case 
can be done in a similar way. 

Similarly to the Section HI the fundamental solutions (p 3 jy(x,uj) are given by [14] : 



<^(x, to) = cos(to(rj — x)) + / G^Tj — x, Tj — s) cos(to(Tj — s))ds, 

J x 

where Gj\f(x,y) are C 1 with respect to the both coordinates. 

The above equation yields to the following representation for the product functions: 

tptf(x, a;)v?j^(x, to) = cos 2 (to(rj — x)) H — / Kfofa — x, Tj — s) cos(2a;(r 7 - — s))ds (25) 

2 J x 

where Ky\f(x,y) is a Volterra kernel, i.e. Kj(x,y) = if y > x and 
Kfox, y) = 2[Gy\ f (x, x-2y) + G^(x, x - 2y)\ + 

/•x px—2y 

+ / Gy\f(x, s)G J j^(x, s — 2y)ds + / Gj\f(x 7 s)G J j^(x, s + 2y)ds. 

Jly—x J —x 

At this point, we note that, as Kjj-(x,y) is a C 1 function, we have 

i J 3 K j H {jj - x, Tj - s) cos(2tu(r i - s))ds = O f^j , x 6 [0, Tj}. (26) 
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Applying the characteristic equation (1251) and developing the products ipjy(x, uS) 
(p 3 j^{x,ui) by (1251) . and finally using (1261) . we have: 

52 ( TTcos 2 (W fe ) ) / ^(s) cos 2 (u;(r i - s))ds = O ( — ), 
o=i Wi J Jo W 



N 



j=i \fc^i 



e (n c ° s2 (^)j \[ (27) 

In the last passage, we applied the fact that J* 3 qj(s) cos2(c<j(Tj — s))ds = O (1/uj), since 
q is in ff 1 (r). 

The left side of ( 1271) is an almost periodic function with respect to ou, in the Bohr's 
sense. Following the same arguments as those of the Theorem [1] we obtain 



3=1 \k^j 



Choosing ujj = vr/2rj and noting that, as the parameters {rj)^ =1 are not entire multiples 
of each other i.e. Ti/rj N: 



,71 Ti- 



cos(ujjTi) = cos( — — ) 7^ for i ^ j. 



2 Tj 



For each value ojj, we have 



] J cos 2 (-cj u.) / qj(s)ds = / qj(s)ds = 0. 

fc#j 'Jo Jo 

and finishes the proof of the Theorem [2J 
5. 5. Proof of Theorem 

In this subsection, we consider two potentials q = an d <? = ®f=iQ.ji satisfying 

the assumptions of the Theorem [5J Assuming that they give rise to the same Neumann 
reflection coefficients, R^{u) = Rtf(u), we have the characteristic equation (T23j) . 

Let us define the operator to be the operator —-£2+ qj(x) over [0, Tj] with the 
domain 

D(C^) = closure of C#(0, Tj) in H 2 (0, Tj ), 

where Cj^(0, Tj) denotes the space of infinitely differentiable functions / defined on [0, Tj] 
satisfying Dirichlet boundary condition at and Neumann boundary condition at Tj. 
By the perturbation theory for linear operators [11] . there exists e\ > small enough 
such that, if |M|z/>°(r), ||<?||l°°(T) < £i then the eigenvalues of £^ remain positive. 
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Considering ((A^) 2 )"^ (A| > 0) the sequence of eigenvalues of Cj^, (1231) implies for 
each j = 1, ■ • ■ , N, 

n^(0>M)^(0,Af) r%(x)<pb{x,)i)<&{x,>i)dx = 0, V< = 1,2,... (28) 

where we have applied the fact that ^^(0, A]) = 0. 

At this point, we will use the assumption Bl' on the lengthes Tj to obtain a Lemma 
on the non-overlapping of the eigenvalues for different branches: 

Lemma 3. Under the assumptions of the Theorem^ there exists e 2 > small enough 
such that, if \\q\\L°°(r), ||<?||L°°(r) < £2, then 

J] ^(0, \i)fa(0, K) ± 0, Vj = 1, • ■ ■ , N, V* = 1, 2, ■ • • . 

Mi 

Proof. In order to prove this Lemma, we only need to show that Xj is not an eigenvalue 
of £^y- nor for k ^ j. 

In this aim, we first show that there exists M > such that for integers ix,t2 > M, 
is different from Af and Af 2 the eigenvalues of and Assume, contrarily, that, 
for all M > 0, there exists z'i,Z2 > M, such that 

A^ = A^ or A£ = A*. (29) 

As the potentials g and g are if 1 , we have the following asymptotic formula's for the 
eigenvalues (see, for instance, Remark 1.1.1, page 7, j5]): 

A , = (2t - i)tt + /j^c^a + c /l 

2r fc 2ir fc \i 2 

- xk _ (2i-l)n | fft(# 1 O (- ) ■ 
2r fe 2ir fe U 2 



However, as by assumption B2, the integrals qj.(s)ds and qk{s)ds are zero, this 
implies: 

Therefore, the relation (129|) implies that, for all M > there exists > M such that 

Assuming, without loss of generality, that ii < i 2 and dividing the above inequality by 
(2ii — l)7r/27fc, we have 



(2ii - 1)tt (2i 2 - 1)tt 
2r 7 - 2r k 
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Therefore, we must have the existence of an infinite number of integer couples i 2 ) 
satisfying the above inequality. However, by Thue-Siegel-Roth Theorem [16], for any 
irrational algebraic number a, and for any 5 > 0, the inequality 



\a-p/q\ < l/\q 



2+6 

j 



has only a finite number of integer solutions p,q (q ^ 0). This, trivially, leads to a 
contradiction and therefore there exists M > such that for i±,i2 > M and j ^ k, 

Mi ± M 2 and Mi * Mr 

For the M first eigenvalues on each branch, we apply the perturbation theory for 
linear operators [TT]. Having j, the branch index, fixed, we will show that for e\ > 
small enough, if H^Hl^, H^IIl^ < e i then the quantities (\DfLi do note coincide with 
the quantities (Af)^ 1; where k ^ j. In fact, for the case of q = q = 0, this is, trivially, 
a consequence of the fact that the branch lengthes Tj and are 2-by-2 Q-linearly 
independent. Now, adding a small perturbation, q or q, in the generalized sense (see 
page 206 [IT]), by the continuity of a finite system of eigenvalues (see page 213 [H]), this 
claim remains valid. Thus, there exists a small enough e 2 = min J= i j ... > jv(e^) such that, 
if lklUoo> ll<?IUoo < e 2 then (A^)f£ 1 do not coincide with (Af)^ , for all j, k — 1, • • • , N 
verifying j ^ k. □ 

Applying Lemma [3] to the Equation (J28|) . for || g|| ioo , || g|| ioo < e := min(ei,e 2 ), we 
have: 



fy(xWjf(x, \l)<ffo(x, %)dx = 0, \/j = l,---,N, Vi = l,2, 
Using the integral representation with a Volterra kernel (1251) . we rewrite 

jT' ^cos 2 Xj (tj -x) + £ K^Tj-x^j-t) 2A ,-^ _ ^ \ ^ = 
^ ' ^-(jc) ^cos2A^'(T i - x) + y ' Tj ' - - cos2A|(r j - = 

So + / ^ 7 " J '~2 <,Tj '~ a?) ^^ d!C J cos2A K r i " *) dt 

where we have applied the assumption B2, qj = 0, for the passage from the second 
to the third line and the Fubini Theorem for the last passage. This implies that 



/ 3 A j (t) cos 2A|(r i - t)dt = Vi = 1, 2, 
Jo 



where 



o 
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By the assumption of Theorem [21 the potentials qj{t) are zero on [0, Tj/2] and 
consequently A^(t) = on [0, Tj/2]. Therefore, we have 

I 3 A 3 {t) cos2\ 3 (r j - t)dt = Vz = 1, 2, • • • . 

J Tj/2 

However, by the Proposition 1.8.6 of [5], the system {cos(2Aj(rj — .))}j=i,2,-- provides a 
Riesz basis over L 2 [rj/2, Tj]. Consequently, 

q.(t) + f '-^ J ' ] qM )<lr = t G [0, 

Jo * 

This homogenous Volterra integral equation has only the trivial solution qj = on [0, Tj]. 
This, trivially, implies that for e = min(ei,e 2 ), if |M|,L°°(r)> ||<?IU°°(r) < e then q = on 

r. 



5.^. Proof of Theorem^ 

We consider two potentials g = ®jLi9j an d <? = ®f =1 qj, satisfying the assumptions 
of the Theorem HI Assuming that they give rise to the same Neumann and Dirichlet 
reflection coefficients, Rj^(uj) = Rjv(uj) and Rx>(u) = Rd(u), we have the characteristic 
equations (1231) and (|2il). 

We define the operator C 3 V exactly as £jy (defined in the previous subsection) with 
Dirichlet boundary conditions at and at Tj . Still, by the perturbation theory for linear 
operators [TT] . there exists e\ > small enough such that, if ||qi|L°°(r), ||?|U°°(r) < e i 
then the eigenvalues of C 3 ^ and C 3 V are all positive. 

Considering ((X 3 ) 2 )Zi (K > 0) the sequence of eigenvalues of jC 3 j^, and ((/x^) 2 )^. 1 
(fi 3 > 0) the sequence of eigenvalues of C 3 V , the equations (T231 and ff2^1) imply: 

J] ^(0, >i)fa(0, A|) P <7,(*)^, Af W M {x, >()dx = 0, 

m .... Vj = l,---,iV,Vi = l,2,..., 

where we have applied the fact that (^(0, A]) = </4>(0> /4) = 0. Following the same 
arguments as those of the Lemma El there exists e 2 > small enough, such that 

k * j 3 = 1, • ■ • , N, % = 1, 2, • • • . 

n^>(o,/4)^,(o,/4)^o, 

Therefore, for each j — 1, • • • , iV we have , 

< &(0>A0> K) >L 2 (0,r J ) = < 9j(0>¥4>(-»^i) >i 2 (0,r,)= 0,Vi = 1, 2, • • • . (31) 
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At this point we define the sequence (v>i)j=i,-,N, t=o,i,2,- as follows 

^ (x) =(0,.-.,0, ^ ,0,..-,0), xGT 
j-th position 

4-i(*) = (o, • • • , 0, ^{x, \i)0ir(x, xi) - i 0, • ■ ■ , 0), x G r, i > 1 

V v ' 

i-th position 

4( x ) = (0, • • • , 0, <f? v (x, Hi)<Pv( x > A) ~ \> °> " " " > °)> x G T, i > 1 

V v ' 

j-th position 

As /J- 7 qj(s)ds = by Assumption B2, (jHTl) implies 

< ?(•), wj(-) >L 2( r)= 0, Vj = 1, ■ • • , N, Vi = 0, 1, 2, • • • . 

where 

N 

<V(-),W(-) > L 2 {T) = ^ < V j (-),W3(-) >L2(0,r J ) • 

i=i 

However, by Lemma [5] proved in the Appendix, the sequence (Wi)j=i,-,jv, i=o,i,2,- forms 
a Riesz basis over L 2 (T) and therefore g = 0, which finishes the proof of Theorem HJ 
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Appendix A. Riesz basis properties 

The goal of this appendix is to provide a Lemma [5] on the Riesz basis property for 
the sequences (t^i)j=i,—,N, i=o,i,2,~ defined in Subsection 15.41 Note that this is a direct 
consequence of a result already proved in [3]. However, we provide a proof for the sake 
of completeness. 

We first provide a classical result on Riesz sequences [SJ page 91]: 

Lemma 4. Let {fj}jL 1 be a Riesz basis for the Banach space B. Let {gj}'jL 1 be 
quadratically close to {fj}f = i, i.e. 

oo 

Then if the sequence {gj}jZ 1 is oo -linearly independent or if it is complete in B, then it 
is, also, a Riesz basis for B. 

Throughout the Appendix, we consider the sequence (ul) to be defined as in 
Subsection 15.41 We therefore have the following lemma: 
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Lemma 5. Assuming q and q in H 1 ^), the sequence (u J i )j = i t ... t N ! j=o,i,2,- provides a 
Riesz basis over the Hilbert space L 2 (T). 

Proof of Lemma{5[ Fixing the branch index j, we prove that {ul(x)}^ is a Riesz basis 
over L 2 (0,Tj). Here, for simplicity sakes, we have identified the vector u{(x) with its 
j-th component and therefore u\(x) is not a vector anymore but rather a function in 
L 2 (0, Tj). We also remove the branch index j. 

First, we prove that {vi(x) := f^-^u^x)}^ is a Riesz basis over L 2 (0,r). Note 

that 



T d t d 

Vi(0) = — -j-(<Pjv.2>)(x, \)<Pm,v(x, \) + — T-(<Ptf,v)(x, \)V>N,T>{x, Aj) 
m dx ' x=o m dx 

t d t d 

nr dx x=r m dx 



= 

x=0 



0. 



In fact, one can easily see that the functions are solutions of 

PX 

- v'l + 2{q{x) + q{x))vi + / N(x, s)vi(s)ds = 4A 2 ^, x E (0,r), (A.1) 



where 

iV(x, s) = g'(a;) + + (q(x) - q{x)) I (g(£) - g(0)d£. 



In particular, we obtain that the functions {vi(x)} c *i l are the eigenfunctions of (IA.ll) 
with boundary conditions Uj(0) = Vi(r) = 0. 

Then there exists a bi-orthonormal sequence {v*{x)}°^ l in L 2 [0,r], eigenfunctions 
of the adjoint operator 

-{v*)" + 2{q{x) + q{x))v* + / N{s,x)v*{s)ds = 4A 2 <, v*(0) = u*(r) = 0. 



Thus are cu-linearly independent. 

From the integral representation of fundamental solutions and from Lemma [U we have 

iirx ( 1 \ 

Vi(x) = sm( ) + O I - 1 as 2 ^ oo. 



So 

oo 



K(ar) -sin(^)|||2 (0jT) < oo 



E 

i=l 

i.e. {vi(x)}'?i 1 is quadratically close to {sin(— By the virtue of the Lemma HI 
this gives that {v^x)}^ is a Riesz basis over L 2 (0,r). 

Now let us show that the sequence of the functions {u^x)}^ is complete in L 2 (0, r). 
Indeed, for any / e L 2 (0,r), suppose that 

f(x)ui(x)dx = 0, i = 0, 1, 2, • • • . 
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In particular, considering it , we have 



f(x)dx = 0, 

Jo 

and therefore, 

/ f(x)(ui(x) + l/2)dx = 0, i = 0, 1,2, 
Integrating by part, we have 

\ X = T 

f(s)ds) (u t (x) + l/2) 



v i( x ) I / f( s )ds ) dx = 0, z = 1, 2, • • • . 



x=0 •'O 

and therefore using the boundary condition ■Uj(O) + 1/2 = 0, 

v i( x ) I / f( s )ds ) dx = 0, i = 1, 2, • • • 



The sequence {"^(x)}^ is complete, hence fj f(s)ds = for x G [0, r]. Therefore 
{■^(x)}^ is complete in L 2 (0,r). Since {r]iUi(x)}°l (where rj^s are appropriate 
normalizing constants) is quadratically close to {cos(^)}°i , it follows from Lemma 
HJthat {ui(x)}°Z is a Riesz basis for L 2 (0,r). □ 
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